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ABSTRACT 


The primary purpose of the Transonic Turbine Test Rig of the United 
τες Naval Postgraduate School, Monterey, California, is to permit 
investigation of turbine performance and flow conditions for transonic 
rotor inlet velocities, This turbine also has provisions for investi- 
gation of the effects on performance of widely varied axial and tip 
clearances, High turbine total-to-total pressure ratios are achievable, 
since an ejector-type exhauster assembly can maintain a partial vacuum 
in an airtight turbine discharge hood, This study is concerned primarily 
with initial turbine performance tests--of the nature of suitability 
tests--both with and without the rotor installed. The installation is 
described, and results of flow rate calibration and labyrinth seal 
leakage tests are included, along with turbine performance test results. 
Also] presented are a simplified three-dimensional turbine performance 
prediction and a one-dimensional analysis of exhauster performance, The 
turbine performance prediction method uses the so-called equivalent 
enthalpy (and corresponding equivalent total thermodynamic properties) 
to permit analysis of the relative flow through the rotor in a manner 


analogous to that used for absolute flows. 
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Symbols 


Fb 


= — 09 


G^" Je ga 


TABLE OF SYMBOLS 
Cross-sectional area (in) 


Constant dependent on square-edged orifice geometry (dimension- 
less). 


Exit opening of a cascade (in.). 

Constant dependent on square-edged orifice geometry (14.71), 
Subscripted, one through three; see Eqs. (122)-(124). 

Blade's departure from being straight-backed (in.).; see Fig. 53. 
Between-teeth chamber heighth, labyrinth seals (in.). 
Constant in square-edged orifice equation ((*R/in.Hg-in.H50)2 lbm/sec,) 
Conversion factor, 2gJ (ft2-1bm/sec?-BTU). 

Conversion factor, 2g]c, ( ft2/sec2-?R), 

Chord of blade (in.). 

Specific heat, constant pressure (BTU/1bm-°R). 

Diameter (in.). 

Kinetic energy (ft-lb). 

Curvature of the back of a blades see Eq. (288) and Fig. 53. 
Force (1b). 

Correction factor, see Eq. (200). 

Universal gravitational constant (32.174 lbm-ft/lb-sec?). 
Total enthalpy (BTU/lbm). 

Energy parameter, boundary layer (dimensionless); see App. I. 
Static enthalpy (BTU/lbm). 

Blade heighth (in.). 

Differential pressure (in. н,0). 

Horsepower, 


Integrand, 
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Symbols 


kis 


= ж ге 


BH. 


Incidence angle (deg. or radians). 

Positive stalling incidence angle (deg. or radians). 

Unit vector. 

Conversion factor (778.16 ft-1lb/BTU). 

Dimension for determining blade-back curvature (in.); see Fig. 53. 


Discharge coefficient, flow metering or labyrinth seal (dimension- 
less). 


Tip clearance (in.). 

Head coefficient (dimensionless); see Eq. ( 81). 

Distance between stations ahead of and after a rotor (in.). 
Mach Number (dimensionless). 

Moment (ft-lb or in-lb), 

Mass flow rate (slugs/sec). 


Exponent used in boundary layer calculations (dimensionless); 
see р, 45. 


Rotational speed (RPM). 

Number of throttlings, labyrinth seals (dimensionless). 
Polytropic exponent (dimensionless); see pp. 79-80. 
Outward normal (unit)vector. 

Pressure (in. Hg. abs). 

Pressure (psia). 

Gas constant for air (53,345 ft-lb/lbm- R). 

Reaction force (1b). 

Reynolds Number (dimensionless). 


Labyrinth seal overall (inlet total to exit static) pressure 
ratio (dimensionless). 


Radius (in.). 


Theoretical degree of reaction (dimensionless); see Eq. (80). 
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Symbols 


te 


tte 


E Ж < с 


Σο 


Entropy (BTU/lbm-?R). 

Non-dimensional entropy (s/c) or entropy increase ( As/c,). 
Cascade blade pitch (in.). 

Surface area (in.?). 

Labyrinth seal between-teeth spacing (in.). 

Temperature (^R). 

Temperature (^F). 

Blade thickness (in.). 

Labyrinth seal tooth thickness (in.). 

Blade trailing edge thickness (in.). 


Projection of Blade trailing edge thickness on the exit plane of 
the blade row (in.). 


Peripheral velocity (ft/sec). 

Velocity within a boundary layer (ft/sec). 
Absolute velocity (ft/sec). 

Relative velocity (ft/sec). 

Flow rate (lbm/hr). 

Flow rate (lbm/sec). 


Fraction of the total flow rate which passes between the hub 
streamline and any other streamline (dimensionless). 


Reference flow rate parameter (in.2); see Eq. (16). 


Non-dimensional radius, r/r,, where ry is the mean streamline 
radius (dimensionless). 


In square-edged orifice equations, Re(1079), (dimensionless). 


In exhauster performance analysis, the axial distance down- 
stream of the inlet to the diverging portion of the exhauster 


piping (in.). 
See Eq. (268), (dimensionless). 


Non-dimensional axial velocity ratio, V,/Va,,, where Уаш is the 
mean streamline axial velocity. 
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Symbols 


п 


In flow measurement equations, expansion factor (dimensionless); 
see Eq, (159). 


Pressure loss coefficient used in determining the loss coefficient, 
Number of blades (dimensionless). 

Absolute viscosity in centipoises. 

з 1007! (centipoises). 


Axial distance in cylindrical coordinates (in.). 


Greek Letters 


οί 


= 


Э, Хул 9 


зо 


ж 
* 


9% оь01037 ОЗ оз ол 


х Ф з = 


Discharge coefficient, labyrinth seals (dimensionless). 

Stator gas flow angle (deg. or radians) 

Coefficient of thermal expansion, flow metering (dimensionless). 
Rotor gas flow angle (deg. or radians). 


Rotor geometric blade angle, 


Ratio of orifice or nozzle minimum diameter to upstream pipe 
diameter, in flow measurement (dimensionless). 


Specific heat ratio, c,/c, (dimensionless). 

Velocity carryover factor, labyrinth seals (dimensionless). 
Labyrinth seal tooth radial clearance (in.). 

Boundary layer thickness (in.). 

Referred pressure Pto/ P. eq» Р-ға = 14.7psia (dimensionless). 
Boundary layer displacement thickness (dimensionless). 
Boundary layer energy thickness 


Less coefficient (dimensionless). 


Correction factor for Reynolds Number, square-edged orifices 
(dimensionless). 


Efficiency (dimensionless). 
Non-dimensional distance from the wall in a boundary layer (y/§ ). 
Referred temperature, Τε Τὰ» T, '47518.4 R (dimensionless). 


Curvature factor (dimensionless); see Eq., (258). 
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Greek Letters 


^ 
^ 


> 


miX 


< нө 5 ντο 


= 


Angle of flow in a meridional plane (deg. or radian). 


Factor used in estimating blading secondary loss coefficients 
(dimensionless): see App. II. 


Pressure loss coefficient, blading leaving loss (dimensionless). 
Area restriction factor (dimensionless); see Eqs. (21) and (187). 
Density (slugs/ft?). 

Non-dimensional flow function; see Eq. (10). 


Labyrinth seal pressure ratio function (dimensionless); see 
Eq. (164). 


Absolute flow velocity coefficient, V/Vpy, (dimensionless). 


Relative flow velocity coefficient, W/Wry, (dimensionless). 


oy Angular velocity (radians/sec). 
Subscripts 

a Ambient conditions, 

a Axial, 

ACT Actual (or computed) value. 

b Refers to flow through stator blade row, exhauster performance 
analysis. 

ο Subsonic compression, in exhauster performance analysis. 

ο In turbine equations, refers to the flow director ring shrouding 
the rotor. 

d In exhauster performance analysis, turbine hood conditions. 

E Refers to so-called equivalent thermodynamic properties. 

e In exhauster performance analysis, supersonic expansion. 

e In turbine equations, entrance station of plenum, 

eff Effective. 

f Due to or corrected for friction. 

H Hub. 
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Subscripts 


I 


i 


15 


ref 


req 


In the vicinity of the cover plate (see Fig. 6). 


Used with efficiency, N, to designate overall turbine (internal) 
efficiency. 


Used with M to denote the moment about the turbine centerline 
acting on the stator-plenum assembly, 


Value for isentropic expansion from inlet conditions. 
Refers to the turbine plenum labyrinth seal. 

In exhauster analysis, downstream of a shock. 

Mean streamline. 

In exhauster analysis, refers to ejector nozzle flow. 
Refers to flow measurement nozzle, 

Denotes inlet conditions. 


Used with the coefficient of thermal expansion to denote the 
square-edged orifice value, 


Used with force, F, to denote pressure forces. 

Plenum condition. 

Used with efficiency, N , to denote polytropic efficiency. 
Used with loss coefficient, $ , to denote profile loss coefficient, 
Relative flow value, 

Rotor value. 

Radial. 

Referred value, 

Required. 

Stator. 

In flow rate equations, refers to shaft labyrinth seal. 


As second subscript denotes value for isentropic expansion from 
equivalent total conditions. 


Used with loss coefficient to denote the secondary loss coef- 
ficient, 
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Ss, 


Subscripts 


T Tip. 
T In exhauster performance analysis, turbine. 
TH Theoretical value, i.e., for isentropic expansion, 
u Tangential., 
u In exhauster performance analysis, upstream of a shock. 
2 Refers to the axial direction, cylindrical coordinates. 
Ө Refers to angle coordinate in cylindrical coordinates 
1 Station between the stator and the rotor, 
l Station at exhauster nozzle exit plane, 
1 Streamline number one (hub). 
2 Station at the rotor exit, 
2 Station at the inlet to the diverging portion of the exhauster, 
2 Streamline number two (between hub and mean streamlines). 
3 Station at the exit of the exhauster piping. 
3 Streamline number three (mean). 
4 Streamline number four (between mean and tip streamlines). 
5 Streamline number five (tip). 
Superscripts 
х Denotes choked-flow, 
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l. Introduction. 

Little qualitative information is available concerning turbine 
performance under conditions of transonic rotor entrance velocities and 
extreme total to static turbine pressure ratios. Additionally, axial and 
radial clearances in turbines for use in rocket applications and gas 
turbine plants are often subject to wide variations because of extreme 
temperature variation, 

Information regarding radial clearance effects on efficiency avail- 
able to designers indicates that tip clearances should be limited to about 
one percent of blade heighth. However, the smaller a turbine, the more 
expensive it becomes to achieve the close tolerances required to maintain 
this percent tip clearance. Also, temperature effects become more critical 
for smaller turbines, with respect to maintaining the desired radial 
clearance over a range of operating conditions. 

Much less information is available to designers on the effects of 
axial clearances on performance. Methods of analysis of turbine perform- 
ance generally used in turbine design assume axisymmetric, uniform flow 
conditions between blade rows, but do not consider the axial distance 
required to permit the flow to achieve such stabilized conditions. 
Consequently, axial clearances between blade rows are dictated primarily 
by space and weight considerations, rather than by knowledgeable consider- 
ations of axial clearance effects, It is even conceivable that optimum 
axial clearances, from a performance viewpoint, may be less than normally 
used for some types of blading. 

The Transonic Turbine Test Rig installed at the Propulsion Laboratory 
of the Department of Aeronautics, U.S. Naval Postgraduate School, Monterey, 


California, was designed by Professor M. H. Vavra of that department to 
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permit investigation of the above-stated problems, This thesis is 

concerned with describing the installation and its capabilities, calibration 
tests, and initial performance tests of the turbine rig, with a stator 
consisting of converging nozzle type blading and an impulse type blading 
rotor, 

Funds for designing and building the Transonic Turbine Test Rig, as 
well as operating funds, were furnished by the Bureau of Naval Weapons 
(RAPP-14). The generosity of that Bureau is greatly appreciated, The 
author is most appreciative, too, of the guidance and counsel provided 
by Professor Vavra, and of the assistanoslin performing experiments so 
willingly given by Mr. L. T. Clark of the Department of Aeronautics. 

2. Installation. 

The Transonic Turbine Test Rig (TTTR) installation is shown schemati- 
cally in Fig. l, whereas Fig. 2 is a photograph of the turbine hood and 
exhauster assembly (with the exhauster rotated 180? from its actual 
configuration in the TTTR installation). The turbine proper is a single- 
stage, impulse-type design. The salient components of the installation 
shown in Fig. l are the: 

1. flow measurement nozzle; 

2. removable standard orifice installation; 

3, plenum labyrinth leak rate measurement standard orifice (two- 
inch pipe orifice, Fig. 1); 

4. turbine plenum labyrinth seals; 

5. shaft labyrinth seals; 

6, turbine hood; and the 

7. turbine proper and turbine plenum. 

This turbine can be operated exhausting either to atmospheric 


pressure or, with the hood and exhauster installed, to a partial vacuum, 
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Supply air is furnished by the Allis Chalmers, VÀ 312 compressor installed 
at the Propulsion Laboratory. Supply air temperature is regulated by 
means of an aftercooler installed downstream of the compressor. Тһе 
maximum temperature of air supplied by the compressor is 300 Ε. However, 
since the turbine has aluminum blading, a limiting value of 200 F has 
been established for turbine inlet total temperature, The maximum 
turbine inlet total pressure attainable is about 43 psia. Exhausting to 
the atmosphere, the total to static pressure ratio possible is about 2.8, 
and a maximum total to total pressure ratio of about 5.6 is anticipated 
when exhausting to a partial vacuum (see Section 4, Exhauster Performance 
Analysis). . Maximum turbine RPM is 20,000 and maximum horsepower is 200. 

The flow rate through the turbine blading during operation is the 
flow measured by the flow nozzle less the plenum labyrinth seal leakage, 
The flow nozzle, with pressure and temperature tap locations indicated, 
is shown in Fig. 3, and the plenum labyrinth seal configuration is shown 
in Fig. 4, Since there is insufficient piping upstream of the flow noz- 
zle to ensure undisturbed flow through the nozzle, calibration of the flow 
nozzle was required. Тһе removable standard orifice assembly is provided 
for that purpose. Plenum labyrinth seal leak rate can not be measured 
during turbine operation; hence, the two-inch pipe orifice is provided 
for determining this flow rate in separate tests. During such tests, all 
flow to the plenum passes through the two-inch pipe and the only exit for 
the flow is through the plenum seals, since the stator exit is blocked by 
,a special cover plate, 

In analyzing exhauster performance, the flow rate to the exhauster 
ducting (see Fig. 5) consists of the flow from the turbine hood and the 
flow through the exhauster nozzle. The former consists of the flows 


through the turbine blades, through the plenum labyrinth seals, and through 
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the shaft labyrinth seals (from the atmosphere to the partial vacuum in 

the turbine hood). There are no provisions for measuring the leakage 

through the shaft seal, and this flow rate must be estimated. Also, flow 
through the exhauster nozzle can not be measured directly. However, with 
the entire flow of the compressor routed to the turbine rig, the flow 
through the exhauster nozzle is the flow rate delivered by the compressor, 
determined as specified by Уауга, (911 minus the flow rate measured by the 
flow measurement nozzle. Additionally, since the flow through the exhauster 
nozzle will normally be choked, the flow rate can be determined quite 
accurately by measuring the total conditions ahead of the nozzle, 

Both the removable and the two-inch pipe orifice assemblies conform 
to specifications for standard square-edged orifice installations as set 
forth by Stearnes, et. al. [2] Detailed information regarding specifi- 
cations applicable to these assemblies, including pertinent dimensional 
information, is given by Eckert. [6] 

The turbine proper and the turbine plenum are shown schematically in 
Fig. 6. The rear support and front hanger shown are attached in turn to 
a cradle flexure assembly affixed to the rig table. The set screw at the 
rear support is provided for setting axial clearances. Figs. 7 and 8 are 
photographs of the turbine with and without the rotor, respectively. Fig. 
9 shows the manner in which the floating armature assembly fits into the 
plenum, As shown in Fig. 6, the following instrumentation is installed 
in the turbine and plenum assembly; 

1. Static pressure (рр) and static temperature (T5) in the plenum; 

2. Total pressure (р), static pressure (pg) and total temperature 


(Ty) ahead of the stator; 


lrigures in brackets refer to bibliography entries of Section 9, 
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3. Static pressures at the hub (Pry) and at the tip (Pip) of the 
stator in the exit plane; and 

4. Reluctance gages for measuring the axial force and the moment 
about the axis, acting on the stator-plenum assembly. 

The air dynamometer used for measuring rotor torque is shown in Fig. 7. 

Spoked-wheel type flexures are installed at both the front and rear 
bearings of the rotor bearing housing for measuring the axial force on 
the rotor. Strain gages are mounted on both sides of each of the three 
webs of these flexures. Pressure taps are installed around the sides and 
in the top of the turbine hood for measuring turbine total discharge pres- 
sure (реа) during operation with the hood and exhauster installed. 

The thrust and torque flexure shown in Fig. 6 was included to permit 
determination of the pressure (рүү) acting on the cover plate concentric 
with and inside the stator ring. The axial force acting on the cover 
plate will be determined by a flexure in the position shown, The pressure 
(ni) can then be determined from this force and the pressure in the 
cavity behind the cover plate (рін): However, the flexure shown has not 
been instrumented, since stator-only tests showed it to be too critical 
in wall thickness and axial groove size, with respect to desired sensiti- 
vity and adequate strength, In the future, this set-up will be replaced 
by a spoked-wheel type flexure, | 

The instrumentation of the stator and plenum assembly was designed 
to permit using the momentum and moment of momentum equations to deter- 
mine stator performance, For accuracy in such performance measurement, 
the pressures at the hub and at the tip of the stator exit, and the pres- 
sure acting on the cover plate must be known, along with the reaction 
force and moment acting on the assembly. Additionally, the average of 


the hub and tip pressures is a better value for stator exit pressure (pi) 
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than either value by itself. During stator tests without the rotor, the 
pressure acting on the cover plate was measured directly by using taps 
installed inside a specially fabricated cylinder placed over the rotor 
bearing housing. However, such pressure could not be determined for the 
turbine tests conducted to date. 

Тһе stator-plenum assembly can be moved axially a distance of 1.2 in. 
to achieve various stator to rotor axial clearances. Moreover, the rotor 
can be moved axially approximately 0.25 in. for increased range and flexi- 
bility in setting axial clearances, 

Fig. lO shows the turbine gas path and pertinent dimensions for maxi- 
mum rotor blade heighth, The basic configuration shown is that for use 
with the converging nozzle stator guide vanes, The dashed lines indicate 
the gas path when the second set of guide vanes, of converging-diverging 
nozzle type, are installed. It is evident from Fig. 10 that a range of 
tip clearances from the minimum practicable to about 0,6 in. are attain- 
able merely by reducing the rotor outer diameter, Additionally, through 
the use of inserts between the wheel outer diameter and the case, various 
gas path flare angles can be achieved, along with varied radial clearances 
and radial overlaps at the tip, obtainable by varying axial clearance, 

The blading profiles of the two available stators are shown in Fig. 
11, апа the profiles of the two available rotors are shown in Fig. 12. 

The area ratio -- exit opening to throat -- of the converging-diverging 
blading is 1.25, corresponding to an isentropic Mach Number of 1.6. Опе 
rotor blading is double-circular-arc profile, while the other is a profile 
with gradually changing curvature, 

The exhauster head (nozzle plus turbine hood flow inlet section) of 
the exhauster assembly shown in Fig. 5 is shown in Fig. 13. Axial 


locations of pressure taps in the ducting are indicated in Fig. 5. At 
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FIGURE 11 
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FIGURE 12 
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each location, four taps are equally spaced circumferentially and commoned 
for obtaining an average reading. Fig. 13 shows also the locations of 
pressure taps along the outer contour of the exhauster head; at each 
location, three taps are equally spaced circumferentially and commoned. 
The location of the total pressure and total temperature tap upstream of 
the nozzle has not yet been determined, Consequently, the location shown 
is approximate only. 

All pressures are measured in inches of mercury against atmospheric 
pressure as reference, using 96 in, manometers with O,l in. graduations. 
Differential pressures across the flow nozzle and orifice taps are 
measured in inches of water, Temperatures are measured using an ice bath 
as reference, 

3. Turbine Performance Prediction. 
General 

The basic method used for turbine performance prediction is that given 

by Vavra. This method is a simplified three-dimensional analysis based 


on the equation of motion and the equation of continuity, with solutions 


A 
<, 


obtained at between-blade-row locations. In contrast to the approach of 


Хал” 
>л 4 


Vavra, the equations of motion and continuity for relative flows, rather 
than for absolute flows, were used in rotor calculations, This procedure 
simplifies computations, since the predicted relative rotor gas exit 

angles can be used directly in the equation of motion, Therefore, itera-, 
tion based on change in total enthalpy across the rotor is not required, 
Also, loss coefficients, estimated using the method of Ainley, [4,5] are 
obtained in the same manner for the rotor as for the stator. Consequently, 
their application must be in a manner analogous to that used for the stator. 


ὄγαντα, M. H., Aero-Thermodynamics and Flow in Turbomachines (New 
York, London: John Wiley and Sons, Inc., 1960) Chap. 16. 
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Proper application is achieved by introduction of the so-called equivalent p. 
enthalpy (defined later). The assumptions made in the performance analysis | 
are set forth below. | 

l. An infinite number of blades is assumed so that downstream 
conditions do not influence conditions upstream, 

2. At the stations of interest, axisymmetric stream surfaces exist, 
so that all derivatives with respect to O, i. e., in the peripheral 
direction, are zero. 

3. Тһе flow is adiabatic and steady so that along any given stream- 
line there is no change of total enthalpy in the stator and no change in 
relative total enthalpy in the rotor. (Note, however, that since cylin- 
drical coordinates are to be used, gradients of total and relative total 
enthalpy can exist in both the radial direction and the axial direction, 
due to streamline curvature, depending on inlet conditions.) 

4. Since the equations are to be solved for between-blade-row 
geht tons | the increase in entropy due to friction within the blade rows 
can be determined using the loss coefficient and the inlet and outlet 
conditions. However, it is necessary to assume that entropy changes 
along the streamlines between blade rows are negligible. 

Development of Formulas 


The basic equations of motion for absolute and relative flows are: 


vH- Vx (VxV) - TVs (1) 
УН = Мх (Ух +22) * T Vs (2) 
Total relative enthalpy, Hp, can be written: 
2 2 2 -- "эргэ" 
H =h+ 3 - =h, t 30 - e И Ақа) 
1 1 
He = һ + Hs. - 4 (3b) 
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The subscripts indicate: l -- ahead of the rotor; 2 -- after the rotor; 
ІН -- the theoretical value; and 2' -- the conditions after the rotor for 


isentropic expansion from Тек =Не/ср Со ру. For adiabatic flow, the 


ча = 
рс нене 


total relative enthalpy remains constant along any given streamline. Hence, 


equivalent enthalpy, defined as: 


| 2 λ 
= Mes . Me 
H, =Ü, + 22: = H. + 5 (4) 
must also remain constant along each streamline, It follows from (3) and 
(4) that:  .- + 
W. ШІ жә СШ 
2 ^ { 
айы: -2- o я е) (5) 


It is evident that the use of equivalent enthalpy permits analysis of the 
relative flow through the rotor in a manner analogous to that used for 
absolute flows, including the use of normal thermodynamic relations. 


Velocity coefficients are defined as: : 


(O = de for absolute flows (6) 
= f elati fl 7 
м ES or r ive flows (7) 

Then, in terms of the loss coefficients: 
жШйы- 1-5 for absolute flows (8) 
хи-1- ї for relative flows (9) 


Fig. l4 shows the relations between thermodynamic properties, velocities 
and loss coefficients in the form of a typical expansion through a turbine. 
The expansion process is shown for three different streamlines in order 
to emphasize the variation of properties which have to be considered in 
the analysis. 

Figs. l5a ~ 15 c show the coordinate system and sign conventions 
used, Station subscripts, zero, one and two refer to ahead of the stator, 
between the stator and rotor, and after the rotor, respectively. 

The equation of continuity is used in its non-dimensional flow func- 


tion -- (-- form. Written in differential form: 
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шн - dAb - dA [21 [(83* - (23 | (10) 
HAE d [1 - - 4АФ = dA} τι) - ч) ] (11) 


for relative flows, Then: 


dA = f Zadr (12) 


Неге, tat is the blade exit opening (see Appendix II) and the factor, Εξ , 








is an area restriction coefficient which accounts for loss of effective 
area due to the boundary layer on both sides of the flow channel, This 
factor is further defined below, Eq. (12) can be made non-dimensional by 
multiplying and dividing by a ng? the product of mean streamline exit 
opening and radius, Hence: 

dA= $ z $&-dXa«m (13) 


Eqs;(10) and (11) in integral form become, respectively: 


IR _ 
ων ο - Om É lm Е е фах (14) 
зь μΕ I α.ἆ η. [ -ᾱ- £ b dX (15) 


( 





Defining: 
| LT ΓΒ < 
We = ыы. - [sz xr ефах (16) 
5 to Stator Xu Stator 


For Wpf» the reference flow rate, in square inches, tat and r, must be in 
inches and Pto in psia.. It follows from (14) and (15) that continuity 


is satisfied, provided that: 


HR Fa dX| = М 
(7 T i i t |2,4 Ju ч- 
Rotor 


Еа. (17) neglects the flow area available because of tip clearance. Тһе 


differential area due to tip clearance is; 
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Letting rg denote the inner radius of the flow direction ring around the 
rotor, the integral form of (15) becomes: 


--— 3 Хе 
wil R. а= | 211 Х 
5 9 |. - и " Xr É Am і f am) 


Clnc 


The Z in the denominator of the integrand is required because the entire 
annular area rather than a between-the-blade area results from integration 
of (19). For small tip clearances, sufficient accuracy should result 

from taking X of the integrand as a constant equal to Ха, Хо, or the 
average of the two, Considering tip clearance area, and assuming that 

the values of ф 3 Ё s Pte» and Тер in the clearance area are approximately 


the same as those for the rotor tip, (17) becomes: 


Ху X. 
Re / Pro d. ЖА ee | X 
{ Tte Й Ти 2-8 043 Ц mn | Te Mis DE dX 


= “ОМ _ 
(аһ 2 Ж 


(20) 


Rotor 


Rotor 

In addition to the requirement that overall continuity be satisfied by 

the stator and rotor flow, it is necessary that the flow rate between each 
777 adjacent streamlines remains constant, 


The term 6 іп (12) is:? 


т QE oes 
ОЁЖ > ου UV | ΝΒ (21) 


Неге, 5, 15 the profile loss coefficient, For the stator, using % їп 

(21) corresponds quite closely to the often-used practice of assuming that 

one-half the total loss occurs to the throat, However, for the rotor, 
ὄγαντα, M. H., Problems of Fluid Mechanics in Radial Turbomachines 


(Rhode-Saint-Genese, Belgium: Von Karman Institute for Fluid Dynamics, 1965) 
VKI Course Note 55b, pp. 646-50. 


4< 


there is no such direct correlation, In (21), H'""* is the energy para- 


meter defined as the energy thickness, p 


» divided by the displacement 
thickness, δι. These thicknesses are given by: 
κκ 


5 = &(1- (179 (75 dn] (22) 


& - &[i-t у дл) (23) 


where: 
4-1 


Xe -H i (24) 


In (22) and (23), n , is the normally-used non-dimensional distance from 





the wall, y/s ; also, m= 1/n, where n is the velocity profile exponent 
in: (w/U) - (y/g )l/n, 

The development of (1) and (2) into usable form is quite lengthy. 
Consequently, this development is included in Appendix I, It is evident 
that these two equations differ only by the term of (2), (Wx20»), which 
represents the Coriolis acceleration, Consequently, the derivation given 
in App. I is for (2). The form of (1) used is obtained from the equation 
developed from (2) by substituting: H for Hj; x, for J, s and U=0. Also 
included in App. I are the derivations of (21) - (23), along with a plot 


of H*** for m = 0.15, The final form of (1) is: 


dUsX) — соза Ды кақ, БЕ Eia 2:6) ) det _ 2 tan% JZ- 


d X 
« 
- θείον e Spes dE — [oriri siwa] d (25) 


and that of (2) is: 


41) --ωφδ]-( Zk $r) (53649 - lton/A se - = - 5178, 


TOU In Pas ο 2 oH C He c 529, аа, alee 
— W^ + W v du -| ЖАШЫН sj dX, (26) 


Qm !2 âm !2 
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Heres 


F< 
! 


- Уа/Уад (ог Wa/Wa..) 
X = r/ry 
k = 4 to 6, streamline curvature factor 


= streamline displacement 


“ОУ 
Э 
| 


p. 
! 


length over which streamline displacement occurs 
s* = s/Cp, non-dimensional entropy 

U = wheel speed, ft/sec. 

Ст = 2gJ = (2)(32.174)(778.16) ft?lbm/sec?BTU, converts H i 

to ft2/sec? from BTU/1bm. 

In the curvature term (+ ἔτι ΤΣ)» the plus sign is for curvature whose 
radius of curvature points in the positive radial direction (see Fig. 15a, 
р. 42). 

Total temperature can be used in (25) and (26) by substituting CT, 


for CjH, where: Сә = 2gJcp. The entropy term in these equations is 


evaluated along the line of constant pressure for the blade exit, The 





required formula 1513 
.2 8,2 
= 2 
* 29701-5) 
Differentiating with respect to the non-dimensional radial coordinate: 
| е Y as 
ОБСЕ: — GhHeose — 
ax 7 dx 4 Y và, (28) 


C, - со5 (1-8) 
Since Y? appears in terms on the right side of (25) and (26), these 


equations must be solved by successive approximations, For a solution to 
be valid, it must satisfy both overall continuity and between-streamline 
continuity, as previously mentioned, Overall continuity is controlled 


^Vavra, M. H., Aero-Thermodynamics and Flow in Turbomachines (New 
York, London: John Wiley and Sons, Inc., 1960) pp. 445-447. 
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primarily by the axial velocity, while between-streamline continuity is 





basically a function of the streamline positions assumed. In integral 


form, (25) and (26) are of the form: 


fa Y: = (14 EE (29) 


where Xj is arbitrary, and ln(C7) is the constant of integration. With 


the boundary condition, Y = 1.0 ас Х « 1,0; 


ο- [ 1dX « e 


Then: 
Án C - [ Lal 
and: 
й X f 
= |, тах – | гах 
Еа. (29) then becomes: 


тҮ" = | тах (зо) 


In solving (25) and (26), the sign convention for the flow angles 
must be carefully observed (i.e., angles are positive in the direction | 
of rotation), In addition, when integrating graphically or numerically 
according to (30), the direction of integration must be kept in mind, 
Integrating toward X (1.0, the summation must be multiplied by -1.0. То 


facilitate solution, (30) can be expanded in an infinite series: 
x 
kÍ IdX 2 m n^ nt х 
m | = а Геи |. 
ү= е -14044-44-4-4(046.2: а-4 ах (31) 


Method of Solution 
Loss coefficients and gas flow angles were determined for both the 


stator and rotor using the method of Ainley. [4, 5] In exception to that 
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method, however, a blockage factor due to trailing edge thickness was 
included in angle determinations, Also, leaving losses were computed 
using the procedure given by Vavra.” Computation of flow angles and loss 
coefficients is given in Appendix II, Loss coefficients were assumed to 
be independent of Mach Number and independent of radius, except for the 
influence of radial variation of incidence angle. Angles were determined 
for Mach Numbers of 0.5 and 1.0. Since only the mean streamline angles 
could be determined for M - l, constant differences between values of 
flow angles for the mean streamline and those for the hub and tip stream- 
lines were assumed for all Mach Numbers. Such differences were taken as 
those for M = 0.5. Smooth curves with inflection points at M = 0.75 were 
drawn between the values of angles at M = 0.5 and M = 1, following Ainley. 
[4,5] Only the converging guide vane, full rotor configuration has been 
analyzed. 

Since the actual curvature of streamlines is quite uncertain at 
best, the curvature was assumed to be zero. The first terms of (25) and 
(26) then reduce to: сов2( ) ds*/dX. When this term is combined with 


the sine squared part of the final terms, (25) and (26) become: 


4(8Ү) > — TN f - Я θεία, i: a с cos dH 


d X, dX, Y, Va ОК, 
2 « 
2 | 99 ор E З 45 (32) 


d Y) "T dó, _ 2 ¿e _ 4Um coy sin, 
d x» li ¿ton qx, X; 58: Y; Way, 


— AU. U cosi 4. |, Gos dH; (- GHecosá | dsf (33) 
We, Х Хо | We | dk 


These two equations, along with the continuity equation in the forms of 


“Ibid., Eq. 5(63), p. 89. 
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(16) and (20), are the major equations to be solved in the analysis. 

Conditions ahead of the stator were assumed as: constant total 
enthalpy; constant entropy; and uniform velocity with Tn tangential 
nor radial component, It follows from the adiabatic flow assumption that 
total enthalpy is constant at the stator exit also, and dH/dX at station 
one is then zero. Similarly, due to the assumption of constant entropy 
at the stator inlet, the only gradient of entropy existing after the 
stator is that due to friction through the blading, as given by (28). As 
a result of the assumed axial-component-only velocity ahead of the stator, 
the loss coefficient for the stator becomes constant at all radii and 
equal to the value for zero incidence angle. 

The analysis was programmed for the CDC 1604 computer, using FORTRAN 
60. The program is included as Appendix III. The procedural steps accom- 
plished by the program, PERFORM, are set forth in succeeding paragraphs 
of this section. 

Five streamline radial positions were used, including the hub and 
tip positions. Inputs for the program are listed below, Except as indicated, 
similar data is read into the computer for both the stator and rotor, 

1, Radii of the five streamlines, in inches. 

2. Values of gas efflux angles in radians for the mean streamline, 
and for Mach Numbers of: 0.5, 0.7, 0.75, 0.8 and 1.0. 

3, Values of the differences in flow angles from the mean stream- 
line value for the hub, mean (zero difference) and tip streamlines. 

4. Rotor loss coefficients and rotor profile loss coefficients for 
five values of the ratio of incidence angle to stalling incidence angle 
(1/15): 0, 0.5, 1.0, -0.8, апа -1.5. 

5. Values of the throat opening dimension for the hub, mean and 


tip streamline radii. 
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6, The number of sets of values given in (7.) below for which 
computations are to be accomplished, 
^ le Conditions for which the analysis is to be done: Mach Number at 
the mean streamline at the stator exit (Mlp); stator inlet total tempera- 


ture (°R) and total pressure (psia.); Mach Number assumed for the mean | 


streamline at the rotor exit; and turbine RPM. 
Considerable use of parabolic equations is made to facilitate compu- 
tation of values required in the analysis. Subroutine PARAB determines 


the required constants, a, b, and c, for the general equation: 


Y =a +bx +¥ex’ (34) 


This subroutine solves three equations in three unknowns using appropriate 
determinants, Gas:angles for the mean streamline at the stator and rotor 
exits are represented by (34) for two ranges of Mach Number: 0.5 to 0.75; 
and 0.75 to 1.0. For a given (input) Mach Number, the mean streamline 
angle can be determined. Using the angle differences between the mean 
streamline and the hub and tip values, the flow angles for the hub and 
the tip streamlines are computed. Then, the parabola for flow angle as 
a function of the non-dimensional radial coordinate, X, is found. 
Constants of Eq. (34) are also determined for: rotor total and rotor 
profile loss coefficients as functions of i/ig, for i/ig - 0 to l and 
i/i = 0 to -1.5; and the ratio of throat dimension, a, to the mean 
streamline value, ag, as a function of X for both the stator and the 
rotor, 

Subroutine STATOR solves (32), without the dH/dX term, by successive 
approximations. Since part of the input is the Mach Number of the flow 
for the mean streamline at the stator exit, and since the flow angle for 


this streamline is known, the axial velocity can be determined from: 
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n= тутун 


lap (35) 
М„= [28 Te (1, - T.) (36) 

and: 
E V, COS e, (37) 


(It should be noted that mass flow rate is yet to be determined.) Іп the 
iterative procedure, the first approximation is obtained with 48 /4Х - 0, 
Іп the second and succeeding approximations, the values of Y determined 
in the preceding iteration are used to find ds*/dX. In determining values 
of Y, using (31), the series through the third power is used. When final 


values of Y are found, conditions at all streamlines are given by: 


2-1! | (38) 
№ = = (39) 
ъ= к – = (40) 
Ιλ. (41) 
P, = Ë Í тат (42) 
V. = MV tane, (43) 


Subroutine FLOWREF is called to solve (16) for the stator, This 
subroutine checks the pressure ratio, р/р» to determine whether or not 
the flow is choked for each streamline. For choked flow, the flow func- 


tion, Ф ; is set equal to $*, where: 


5l 








8-1 et T (44) 


This equation is obtained from the expression for ф in (10), with: 
T 
gem m WS, 
Ro Се ) (45) 


Values of H'** аге computed using the expression derived in App. I: 




















| | Xe xi xe à Xe 

... Хе-| Ч Jm tl + Smtl % 7m +! ч 910141 Ç Itm+! 
| | Хе Хе > Хе | 

Хе-| m+! Jm+l 5m *l 7m*l Әт! 


The area reduction coefficient, f , is then found from (21). Integration 
is accomplished, and between-streamline flow fractions are determined. 
With wr defined as the fraction of the total flow which passes between the 
hub and any other given streamline, there is: 
С 
70 1:4Х 


É (47) 
f Iu dX 


ч 
Here, I. is the value of the integrand of (16). The actual flow rate in 
lbm/sec is given by: 

й о шилийн а (46) 

It would be desirable to have one-quarter of the total flow pass 

between each pair of adjacent streamlines. However, this condition is 
not absolutely necessary, since for this analysis zero curvature was 
assumed and conditions ahead of the stator are completely uniform. As a 
result, the conditions at the stator exit can be used to establish both 
overall and between-streamline continuity requirements. Slight differences 


from multiples of one-quarter for the flow fractions given by (47) are 
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then tolerable. It remains mandatory, though, that the between-streamline 
flow be the same for the assumed streamlines at the rotor exit as for the 
corresponding streamlines at the stator exit, 


Subroutine ROTOR] is used to convert the absolute flow conditions 





ahead of the rotor to relative conditions. For angles and velocities: 


Ше Малы — ; fin inches | (49) 
Wa, = Vu, - U, (50) 
Á= tori C (51) 

C (52) 
U; = U -7 (53) 


In (53), ro is the radius of the given streamline at the rotor exit. 


Other necessary values are obtained from: 


_ Wi. ui - U? 
Н; = C, Ti + DN + JN (54) 
he = да, (55) 
апа: 
" 


uae (56) 


The values of Hg, Ttg» and Pc, are for the fluid entering the rotor along 
streamlines chosen for the stator and leaving the rotor along the same 
streamlines at radii assumed for the rotor, Similarly, the values of 
entropy, from which entropy increases through the rotor are determined 
for the several streamlines, are those for the matching streamlines at 
the stator exit, These entropy values are the increases in entropy 


through the stator blading, as measured from the constant (base) value of 
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entropy existent at the stator inlet, Hence, the rotor inlet value of 
entropy can be considered known for each streamline and the gradient of 
entropy with respect to X at the rotor inlet can be found. Also, since 
Hg is constant along each streamline through the rotor, the gradient, 
dH;/dX9, is known, It is convenient to write: 

gb з о 3 дан (57) 


мһеге sjo and ds} ο/άΧ» are known, Then: 





45 wm "EM 
JX apr ЧО er) 


The unknown gradient in (58), 421/4хҙ, is found using (28), written in 
terms of the relative flow properties. 

The rotor loss coefficient required for the solution of (28) is 
determined in ROTOR] by using the parabolic functions determined earlier 


in the program, The independent variable, i/i,, is given by: 


3 E (59) 


Here, ES. is the blade angle. The resulting value of i is divided by 

the positive stalling incidence angle: (1354 /5753) /ғадіап. Loss coeffi- 
cients, both total and profile, for the rotor are taken as constant and 
equal to the values for i/i, - l and i/i, - -1.5 should the value of i/i, 
be greater than 1.0 or less than -1.5 respectively. 

Subroutine ROTOR2 solves (33) by successive approximations using the 
values computed in ROTORI and assumed values κο S the mean stream- 
line and of Yo for the other streamlines, Тһе assumed values of Y, for 
the first iteration are estimates. Then, for succeeding iterations, the 
preceding iteration?s values of Y, are used, For convenience, the non- 
variant part of ds] /dX,, is included with the invariant terms: 

-2tan,8, dés - 2 sing. Since the gradient of Y,at the rotor exit is quite 


š 
wy yup MS + э . . DN 
.. %.% % š ς᾽: “> ΚΠ. 


54 





severe, terms of (31) through the fifth power term are used. 

When either satisfactory convergence of Yj values is obtained, or 
the prescribed number of iterations--13--have been complete Таса буынын 
conditions at the rotor exit are determined using (38) - (43) written for 
equivalent properties. The absolute tangential velocity, Үдэ» is obtained 
from (50), written for the values at station two, Satisfactory convergence 
of Y, values is defined as: no change greater than 0,005 in successive 
iteration values of Y, for any streamline, 

Overall continuity is next checked using subroutine FLOWREF with the 
rotor exit properties, Eq. (19) as used in the computer program is writ- 
ten to include the flow through the tip clearance area as part of the 
flow between streamlines four and five. Written for the trapezoidal 


method of integration, the integral of (20) is: 


ο ο πα πω. 


When 1 = 4, the first term in the brackets becomes: 


Βε/β. a 20 6 (A/e.) 2) 61 
2 2 Я17 "(ей Є Am сно mn (61) 


Неге, К is the tip clearance and, therefore, k/r, is the AX, of the tip 





clearance, The term (X5-X4) in the denominator cancels in the solution 
of (60). Eq. (60) must yield the value of ҚАРАСЫ Zr Rotor given by “rf 
determined in the stator solution, If the rotor flow is less than 
required (i.e., the value from (60) is too small), the assumed value of 
axial velocity for the mean streamline is increased and, if greater, the 
mean streamline axial velocity is reduced,  ROTOR2 and FLOWREF are then 
called successively until continuity is satisfied. 


When overall continuity is satisfied, subroutine SLINE is called to 
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check between-streamline continuity. This subroutine checks the fractions 
of flow between streamlines one and two, one and three, and one and four 
at the rotor exit against the values of flow fractions between the same 
sets of streamlines at the stator exit. If the rotor and stator fractions 
for all three sets are equal within tolerance, the solution desired is 
known. If inequalities exist, the radial locations of streamlines two 
three and four at the rotor exit are adjusted and successive approximations 
continued, Also, all variables required in ROTOR2 and for the rotor 
equations in FLOWREF are adjusted, as required, to values appropriate for 
the new streamline locations, Then, ROTOR2, FLOWREF and SLINE are called 
as before until the equation of motion and both continuity requirements 
are satisfied. 

It is appropriate at this point to discuss the methods used to 
determine required gradients, to correct radial locations of streamlines, 
and to correct gradients for new values of Х,, For variables for which 
parabolic equations were determined as f(X) initially, the gradients are 
given by: 
dt). b +2cK (62) 
The parabola itself and (62) need be corrected only for changes in the 
radius of the mean streamline, since the value of the variable at a given 


radius does not change. Then, since; 


X= Xie =f (a) (63) 


ín Іп vi 
where the primes indicate the new values, the resulting correction factors 


required to permit using X? in the equations are: 


— m 
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(64) 
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Бал (65) 
The constants, b and c, of the parabolic expressions of the type (34), 
and occurring in (62), are then: b? = Fob; and c! 7 Fo,c. The constant, 
a, is unchanged. 
New values of streamline values of X = X', in terms of the old value 


of ry are obtained using flow fractions (see (47)) and: 


dX 
TA мен 


Here, the positive sign corresponds to a required increase in X due to 


! у : 
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too low a fraction of the total flow passing between the hub streamline 
and the streamline in question; the minus sign, then, corresponds to too 
great a flow fraction and a required decrease in X. The situation is 
illustrated in Fig. 168. Note that the absolute value of flow fraction 
difference is used in (66) and that the gradient, dX/dwg, is always taken 


for increasing X. Then, in the nomenclature of Fig. l6a, there is: 


ο ο. y XD 
de, 4()-6(0) 121 150) 


The value, X*, obtained yields the new radius to be used, with: 
| 
T. (68) 


Since the hub streamline flow fraction is zero and the tip value is 1.0, 
only the radii for streamlines two, three and four need be corrected. New 


values of X are given by: 


Fa L (69) 


new | ine w) 


New values of 0, are obtained for the new values of гу, using? 


new 
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The curve of Fig. l6b is ыл to represent some arbitrarily 
selected f(X) for illustration of the method used to determine gradients 
of entropy and equivalent enthalpy. Gradients at the hub and the tip are 
taken as the straight-line slopes between values of f(X) for streamlines 
one and two and streamlines four and five respectively. Gradients at inter- 
mediate streamlines (two, three and four) are taken as the average of the 
straight line slopes between values for these streamlines and the values 
for streamlines on either side. All slopes are computed for increasing X, 

The manner in which values of s*, Нр and the gradients of these two 
properties are determined for corrected streamlines is illustrated in 
Figs. l7a and l7b. The curves of Fig. 17a show the possible gradient 
variations for curves which might pass through the property value for the 
streamline of interest, Also shown are the straight line slopes to the 
property values for the streamlines to the left and to the right, Fig. 17b 
shows the increasing X portion of Fig. 17a with curve segments, straight 
line slopes and actual slopes at the streamline value of the property. 

To obtain the rate of change of the gradient, it is assumed that the 
straight line slope is the slope of the curve at the value of X one-half 
way to the adjacent streamline. The rate of change of slope is then the 
straight line slope minus the slope of the curve at the streamline point 
on the curve, divided by one-half the difference in X values between іт 
two streamlines, Note that the positive case is used in obtaining the 
new value of X, whether this value is greater or smaller than the old 
value, If X must be decreased, the straight line slopes to the left of 
the streamline of interest are considered to extend to the right of the 
streamline point, Then, once the rate of change has been determined, the 


gradient for the new value of X for the streamline is the old gradient 
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METHOD FOR DETERMINING RATE OF CHANGE OF GRADIENTS 
FIGURE 17 


plus ( AX є Хпеу 7 Хола) times the rate of change of slope with X. It is 
evident that the correct gradient results, since the sign of AX deter- 
mines whether the plot of Fig. l7b or its mirror image, with respect to 
the old streamline value of X,is used. Similarly, using this new slope, 
the value of the variable at the new X is obtained from the old value 
plus AX times the new value of the gradient. : 
Once the required solution of the analysis is known, additional pro- 


perties and quantities are determined for each streamline, The change in 


total enthalpy is given by: 


AH -H,-H, 7 (U7- UV.) στ. (71) 
Then, sequentially, the following are determined: 

= AH 

el, cw (72) 
шан 

№ = p Й - 
xr. 

«єю (74 


The temperatures, Toig: and Το» for isentropic expansion to p, from Pto 


and Ре respectively, are given by: 


ы ια} (75) 


үз 
Th, = Кт і (76) 


The overall turbine efficiency (internal efficiency) is then: 


С ЛЫШ T. t, 
ni =h. NUS (77) 


Additionally, the stator efficiency, Ng» and the rotor efficiency, 7g, 
are computed using (78) and (79) below, even though they are related to 


the input loss coefficients by (1-8). The values obtained in the output 
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provide a check on program operation, 


x m | 
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Stator and rotor velocity coefficients, Q and (p, are obtained as the 


square roots of (78) and (79), respectively. The theoretical degree of 





reaction, r*, and the so-called head coefficient, kig, are given by: 





Жы This — Thi 
maa A m жү 


ο 215 


4 --- Ah; = 2 cs (Ts —Tzis) qJ 
Κι EE ρα £ Т (81) 


Γ (80) 


In order to obtain mass-flow-weighted values of moment and horsepower, 
the following equations are used, where we is the flow fraction at the 


rotor exit, and the subscript w indicates the mass-flow-weighted value: 


: С te) M) 
(AH), = p Є * AH;) йы ы (82) 
Then: 
— (AH), J ¥ 
К Е (83) 
Mg, = Вэ 50) (84) 


Referred values are obtained, following NASA practice. For Y= 1.4: 





Hre = T (55) 
Meret = 2 (86) 
Nee = s (87) 
175 ш | (88) 
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9505: 70582 (89) 
= Ro 22 Pee 
о = Pro 14.7 (90) 


Results and Discussion 
The performance analysis was accomplished for assumed stator exit 
mean streamline Mach Numbers, M 1? between 0.9 and 1.4 in steps of 0.05. 
The general range of RPM used was 11,000 to 20,000. However, at assumed 
Mach Numbers of l.l and above at low RPM, the high values of rotor loss 
coefficients caused divergence in the iterations for Y, for the tip stream- 
line, Also, for Mn = 1.4 and 18,000 RPM, and for М а - 1.35 and 18,500 


RPM, the pressure ratio across the rotor, / ра» indicated choked flow 


Peg 
іп the rotor and overall continuity could not be satisfied. Since no 
provisions for this contingency are included in the program, 17,500 and 
18,000 RPM were the upper limits for Min = 1.4 and 1.35, respectivelyd 
Table I indicates the values of Mn and RPM for which the analysis was 
attempted without success. 

TABLE I 


PERFORMANCE ANALYSIS INPUT VALUES FOR 
WHICH RESULTS ARE NOT AVAILABLE 
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Stator inlet conditions--total temperature and total pressure--were 
chosen to give a static pressure of 14.7psia for the mean streamline at 
the stator exit, except that limiting values of 660°R and 43psia were 
used, At the higher Mach Numbers, Mj, more than 1.25, stator exit plane 
pressure was a partial vacuum, Fig. 18 shows the schedule of input 
variables as functions of Mj,. 

Representative output from the program is included in Appendix III. 
Since most of the printout column headings are self explanatory, no 
general listing is givens however, App. III includes such description of 
print-outs as is required for clarity. 

Table II summarizes the results of the performance analysis shown in 
Figs. 19 through 30 for ease of reference, 

The analysis showed values of Хү to be essentially invariant with 
М а Consequently, a single curve of Y] versus Xj is representative for 
all values of Μι. Moreover, since the flow discharge angle, o , was 
assumed to be constant for Mj = l, and since &, values for Mj - 0.9 
differ from the M] = 1 values only by 0.05 , a single curve of Vy, versus 
Хі із well representative for all values of Mj). 

Rotor exit conditions vary with each of the rotor inlet condition 
variables, Hence, the curves of Figs. 24 - 26 are included to represent 
the general nature of rotor exit conditions, as given by the analysis. 
However, the analysis showed that rotor exit plane pressure, po, is 
essentially constant with radius for each set of rotor inlet conditions. 
The maximum difference in values of po for the several streamlines was 
approximately 0.5 psia for given rotor inlet conditions. It is evident 
from Figs. 24 and 25 that the method used for determining gradients of 
entropy and equivalent total enthalpy is satisfactory, since there are no 


x ама š 
marked inflections in the curves of s and Hg. It is interesting to note 
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TABLE II 


SUMMARY OF PERFORMANCE ANALYSIS 
CONTENTS OF FIGS. 19-30 
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that the band of variation of Y with RPM for constant Mach Number, Mg 
is relatively narrow. This observation, although valid in general, is 
offered with qualification, since, as stated previously, results could 
not be obtained for certain assumed inlet conditions for which the 
iteration for Y, became divergent. 

Fig. 19 indicates that wWper decreases with increasing pressure 
ratio, Pto/ P1; whereas for choked flow it could be anticipated that Wref 
would be constant. However, it is evident from Fig. 20 that the reason 
for the decrease is that values of the area restriction factor decrease 
as peo/p] increases. 

Considering that normally the rotor exit plane pressure is known, 
selected performance parameters were plotted versus Pto/P2; so that by 
cross-plotting, these performance parameters could be shown as functions 
of N/([8 and k;, for different values of Pto/P2 = constant, The resulting 
curves, presented in Figs. 27 - 30, are designed to permit performance 
estimations for selected values of RPM, inlet total temperature, pressure 
ratio Ρτο/Ρ» and either pc, or pj. Figs. 27 - 30, in conjunction with 
the stator and rotor gas angle information included in App. II, and Figs. 
19 and 21 - 23, provide sufficient information for determination of mean 
streamline performance. Recommended values of mean streamline radii for 
use in such a performance estimation are 4.26 in. for the stator and 4.45 
in, for the rotor, Variation of conditions with radius in the stator 
exit plane can be found using Figs. 20 - 23, Moreover, since ро 15 
essentially constant with radius, an approximation to variation of 
properties in the rotor exit plane can be obtained, by using po = constant 
and estimated values of Y taken from Fig. 26. 

It might appear, on the basis of the curves of Fig. 28, that the mean 


streamline values of internal efficiency are not adequately representative 
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of turbine performance, since the mean streamline curve falls appreciably 
closef|to the tip curve than to the hub curve, However, the average of 
the 13168 of N; Гог all streamlines is generally only 0.001 to 0.002 
different from the mean streamline value. It is significant, too, that 
the values of AH for all streamlines vary little from the value for the 
mean streamline, Consequently, the mean streamline values of internal 
efficiency, as well as of the other performance parameters, can be 
considered excellent indicators of turbine performance as predicted by 
the performance analysis. | 
Unfortunately, time available did not permit thorough investigation 
of the circumstances leading to the indications of rotor choking. Outputs 
from the computer for indicated choked flow in the rotor show that the 
manner in which the area restriction factor and the flow function vary 
with pressure ratio might cause too-early indications of rotor choking. 
As the pressure ratio approaches the critical value, the rate of increase 
for Ф decreases while the rate of decrease of $ is almost constant, so 
that thein product: dnoreases only slightly. It appears that the rate of ^: 
increase of this product becomes zero at the critical pressure ratio and 
then decreases with increased pressure ratio, pep/p2 in this instance, 
It will be recognized that this decrease of the product, ф 6 , was shown 
to account for the non-constant referred flow rate, It is conceivable 
that, should rotor flow calculations show deficient flow rate at pressure 
ratios reasonably below the critical value, any increase in axial velocity 
yields increased pressure ratio with no increase in flow rate, Hence, 
large increases in axial velocity and pressure ratio result with succeeding 
iterations. The propriety of the effect of (D in the foregoing is beyond 
question, However, whether the apparently overriding influence of the 


area restriction factor, Ë , is correct is considered problematical, It 
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is pertinent, though, that for assumed inlet conditions of Mn - 1.35 and 
18,500 RPM, for which the analysis showed the rotor to be choked, if § 
were held constant after the first determination of rotor flow rate, the 
maximum rotor flow rate given by using %” would not have increased the 
rotor flow rate to the value required. In this instance, then, the choking 
of the rotor was not caused by the decrease in values of & . The same 
observation holds also for the other choked-rotor result of the analysis. 
The divergence in Yj values which prevented solutions at low RPM was 
probably caused in part by the higher loss coefficients for the rotor due 
to incidence angles approaching positive stall incidence. However, the 
loss coefficients by themselves should not have strong enough influence 
to generate the divergence. Rather, it is more probable that the initially 
assumed values of Y were enough in error to start a divergence pattern 
after the first iteration, since either Y? jj/Y5 .. or that ratio squared 
is a multiplier of the several terms of the integrand for determining Y, 
(see (33), p. ^48). Considering the effect of Ё іп (28), р. 46, it is 
apparent that the higher the value of 5 s the greater is the subtractive 
term in the denominator. Still, once the value of f has been introduced 
into that term, later fluctuations or changes in the term are due to YZ 
(or to the ratio given above squared as a multiplier), for axial velocity 
considered as a constant, It is then apparent that high rotor loss coef- 
ficients can нь said to establish a certain level of bias in the entropy 
terms which, when aggravated by relatively large fluctuations in Y, values, 
can ultimately drive the subtractive term of (28) to values greater than 
one, Of course, this latter event is the final manifestation of the 
divergence of Yj values, since the logarithm of a negative number is 
undefined. Obviously, the subtractive term of the numerator of (28) is 


affected by Yo value fluctuations, but only about one-half as much, In 
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the instant analysis, the large gas efflux angles aggravate the 'bias'! 
produced by the loss coefficients, since the cos? g in the denominator 
serves to increase the subtractive term, 

The method used for this performance analysis appears sound and is 
deemed relatively easy to use. Convergence of iterations is quite rapid 
except as noted above, It was found that the rate of convergence was 
increased if initial input RPM, for otherwise constant input parameters, 
was the middle value of the expected RPM range, with subsequent RPM inputs 
increased or decreased from the initially assumed value. Obviously, this 
method, like all methods, can be considered no more accurate than the 
assumptions used in the solution. Undoubtedly, the most critical values 
assumed are those for the loss coefficients. 

The program PERFORM used in this analysis, although apparently giving 
satisfactory results, can surely be improved. It was noted that values 
of a/a, for the mean radius slipped somewhat from the desired value of 
1 (to 0.9943 in one instance). Since this ratio is expressed as a para- 
bolic function of X, and since slippage was due to the correction factors 
for the constants in the parabolic expression, the method for correcting 
rotor exit streamline radial locations should be reexamined with a view 
to decreasing the noted discrepancy, Also, it is considered that a more 
appropriate tolerance for flow fractions would be a percentage, say one 
percent, of the required flow fraction, This change would yield more 
uniform requirements for between-streamline flow rates, rather than the 
presently-used two, one, and three-fourths percent tolerances, used in 


order, for the one-fourth, one-half and three-fourths total flow points. 


78 





4, Exhauster Performance Analysis. 


General 

Analysis of the turbine during exhauster operation was ew to 
obtain preliminary estimates of total-to-total pressure ratios achievable 
across the turbine and to predict shock conditions which may occur in the 
exhauster piping. The analysis conducted was for one-dimensional flow, 
based on the momentum, continuity and energy equations. Extensive use 
was made of polytropic process relations and the non-dimensional flow 
function, ф » Which is obtained by combining the energy and continuity 
equations, The relations for expansion and compression processes are 


given below without proof, 9 For an expansion process: 








з = = ae E. τ (91) 
me 44 > ΥΠ m В, 
where: 
nz ART (92) 

01577 polytropic efficiency 

Y - specific heat ratio 

м - flow rate (lbm/sec) 

Tto- stagnation temperature at the inlet (°R) 

A - flow area at the station of interest (in2) 

Ptg- stagnation pressure at the inlet (psia) 

p - static pressure at the station of interest (psia) 

R - gas constant for air (ft-lb/lbm-?R) 

g - universal gravitational constant (1bm-ft/lb-sec2) 


Similarly, for a compression process, there is: 


буауга, M. H., Problems of Fluid Mechanics in Radial Turbomachines 
(Rhode-Saint-Genese, Belgium, Von Karman Institute for Fluid Dynamics, 1965) 
VKI Course Note 55a, pp. 22-35. 
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i Ла a la (a VE (o VR 

i EC 93 
АЧ ies = ф = Pt" ре 
Here, the mn one and two refer to inlet and station of interest 


respecttvély. Тһе polytropic exponent is given by: 


„_ 
ЕЕ πα (94) 


Іп (93), the quantity T is the static temperature and the other variables 
are as defined and have dimensions as specified for the expansion process. 

Fig. 5, p. 26, shows the complete exhauster and its dimensions and 
Fig. 13, p. 37, shows the exhauster nozzle and the flow channel from the 
turbine hood. Fig. 5 also shows station designations and subscripts used 
in the analysis and the general locations where the several polytropic 
efficiencies apply. 

Method of Solution 

Choked flow is assumed for the flow through the stator blade row and 

the flow through the exhauster nozzle. For this condition, the pressure 


ratio in (91) is: 


e. 


Жемісті jo (95) 


Then, (91) becomes: 


E 
фо БА (0 " 


For choked conditions, flow rates must be determined using this value of 
ф- ф*; however, (91) remains valid for determining pressure ratios lower 
than the critical value. 

Fig. 3l shows the processes involved in the analysis as they might 
appear on a temperature-entropy diagram. Turbine and nozzle flow processes 


and the several solutions possible downstream of the nozzle exit plane are 
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indicated, Station designations used in Fig. 3l correspond to those given 
in Fig. 5. At station two, the turbine and nozzle flows are assumed to 

be mixed and the combined flow is assumed to have uniform properties 
perpendicular to the flow direction, 

As previously stated, the flow at station one consists of the flow 
from the turbine hood plus the flow through the exhauster nozzle. Тһе 
former flow rate is the sum of the flows through the turbine, through the 
plenum labyrinth seals, and through the shaft labyrinth seal, Writing 
the sum of these flows, wp, in the order named: 


The equations for Ww, and Ws, as derived in Section 5, are: 


—r —À — — -———— 





А А, 8, 9 ' 3. pe "t / Re 

= § Fg i $= | oni = 
A, P. [4 RE Si 2 

. T 5 2. 5 ο - 7 о. 

ik i ὦ. 22222) - 


The flow rate through the turbine is obtained from (91), where б-ф” 15 


$* оз А 8.9 (100) 


For the expansion process from the hood to the exit plane of the nozzle, 


given by (96): 


it 1s convenient to writes 


ων ο 
E - Ф = E -! | 9 | "d Us Я emp) 


Then, putting (98), (99) and (100) in terms of O's, there is, with (97): 


аА ФА фф А [E (102) 

















Implicit in this relation is the assumption that the total conditions at 
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the inlet to the plenum labyrinth seals are the same as those ahead of 


the stator of the turbine, Combining the left side of (102) with (101): 


p ке wu. нь A, |Ë Р” р АЛЬ! қ ug (103) 





= |+. =Q, A Р СУ; TRIE = ф, А, Тела (104) 


Ра / Ріт 


Here, the exponent, Np, is the value corresponding to the polytropic 
efficiency assumed across the turbine. This expression can be written, 


with (102): 


INO 


мн тан: OA md ns | 2 = GA+OA,+ OA, ЯГ 


т („/в,) ἔπε 


Then, introducing the right side of (101): 








(105) 





np 


PNE lm ын ах |). ыу Е. NN 
(n, /0,) т sal Б. | = fta ын À, + QA, + АЛГИ T (106) 


The form of this equation used in the analysis is obtained, after some 


rearrangement, as: 


7 
E n4 en] Leal 
Lar (aj [ijt πμ fan 


In (107), there remain two unknowns: pj, and Ptg’ Hence, it is 





necessary to solve the equation by iteration, using desired values of Pto 
and Tto’ Although two solutions of (107) exist, the supersonic solution 
is not of interest in the analysis and only the subsonic КЕЕ is 
obtained. With an assumed value of Pip? Pty is determined from (107), 
the condition and velocity of the fluid flowing to the exhauster from 


the hood are established from the relations: 


83 


Пт-і 


nr 
= V UR) Vise» 
: nazi 
Ur n 
The = М in 4 (109) 


VETT OMEN (110) 


Eq. (96) with (91) givesthe flow rate through the nozzle; with the 


asterisk denoting M - 1,0 at the throat: 


| = = 
‚ SAÀPSÍS (2 ще nig. Lo T 
к=к тн (ду [к] чи 


Since the flow rate at the nozzle exit is the same as at the throat, for 








the exit, station ly; there is; 


міт - b RES A _ s 
Αν EE | QUA, т “η 


Here, the exponent ny is used to indicate that different polytropic 











(112) 





AQ | efficiencies are assumed for flow to the throat and from the throat to 

Y | the exit, Eq. (112) can be solved by iteration to find Рік" Although 
both a supersonic and subsonic solution of (112) exist, the subsonic 
solution was not considered, since the back pressure required for a shock 
at the exit is about 27 psia for an isentropic nozzle. With Ply deter- 
mined, (109) and (110), written for the nozzle variables are solved to 
determine the temperature and velocity at station ly. 


The momentum equation, written between stations one and two, 15: 


м d цан хэй, A г s AC P,- Pa) 258 (218) 


The friction force, Fr, can be expressed as a pressure 1055: 
2 
rep Age (114) 


This pressure drop, in turn, can be written as? 
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ар 0-0-0 (115) 


Here, L is the duct length between the two stations, D> is the duct 
diameter, V and jo are the average of the values of velocity and density 
at stations l, and two, and f is the friction factor. Conditions for 
the flow through the outer annulus at station one are used, because the 
pressure loss is considered to be caused primarily by the friction along 
the duct walls. The factor, f, depends on Reynolds Number and pipe 
roughness. The approximate range of Reynolds Numbers obtained from 
preliminary calculations was 1.7 to 2.2 (10°); and, for a surface rough- 
ness of 0.00015 ft. for commercial steel pipe, a surface roughness to 
diameter ratio of 0.00024 was obtained. With these values, f-= 0.017 was 
determined from the Moody Diagram given by Eschbach. [1] 


Letting рі, » p] - AP, (115) becomes: 


W, +W W уу = 
(сег) У, ο η ний А. (Р, - P) (116) 
The total temperature of the fluid stream after mixing is given by: 


, . Г. | | | 
L, = A s £ W Tha š W = Ур (117) 


Then, from continuity, there is: 





à 
ἪΝ, wk (* x | 118 
үен RÅA RA % 2919, а 


Solving (118) for p, and, with (117): 








S WR [wtih . α | 
ae Far (a тәг, ши. 


Next, combining (119) and (116) to eliminate po: 


М " " А WR Mk ГМ - 
Mie — ey AMT (roue) 029 
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Multiplying (120) by Və, expanding the last term, combining like terms 


and transposing, the form used is obtained: 








V. E ET A - V s "Ма * Vi, +A, P, | Tile = O (121) 

Let: 
B= 9 (!~ 255) (122) 
B, =- Hv, + шү, AA (123) 

and: 
В; = МКТ, (124) 


The solution of (121) is then: 


T En BE p 
сийн 17 2-8 - 


This equation will yield both a subsonic and a supersonic solution. With 
the velocity, V5, known, the fluid properties at station two are obtained 


from (119) and: 


ТА 

Fn "M (00 
Th - 

Жа 5») (127) 


For the subsonic solution of (125), the fluid is compressed in the 


diverging portion of the exhauster and (93) applies as: 


42 (ЕБ E Πεν, 
A ШЕ $ = ο μα ў =) 


This equation is solved by iteration to determine рз, Since the flow must 





adjust to atmospheric pressure at the exit plane of the diverging portion, 
if the value of pz obtained does not equal atmospheric pressure, the value 


assumed for Pie, is in error. For p,4pg, a higher value of Pip is assumed 
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and for p3> Pa, a lower value is taken, Iteration is then continued 
until the value of Ply for which p3 = pg is found. Values of fluid 


properties and the velocity at station three are then determined from: 


n. -! 


r=” (129) 

з 2/ 24 То (T, - 13) (130) 
FS 

Pe, = R RIT (131) 


«αν 


The supersonic solution of (125) gives the velocity at station {νου 





| The temperature is obtained from (126), the pressure from (119), and the 

| total pressure from (127). Knowing these properties, the location of a 
normal shock can be estimated. The possibilities which must be considered 
in the one-dimensional analysis are: 

l. a normal shock at the inlet to the diverging portion, followed 
by subsonic compression; 

2. supersonic expansion from the inlet to the diverging portion to 
some axial distance, X, downstream, with a normal shock at that point and 
subsequent subsonic compression; and 

3. supersonic expansion to the exit plane of the diverging portion, 
followed by an oblique shock outside the exhauster, 

Not listed in the above are supersonic expansion to an exit pressure 
either equal to or less than the full expansion pressure. In the analysis, 
this consideration is not of practical interest; however, consideration 

of these possibilities is included in the solution procedure for complete- 
ness, 

Since the flow must adjust to ambient pressure, possibility (1.) is 
checked first, because it provides maximum pressure recovery, ,If the 


resulting value of p3 is less than pg, there can be no normal shock from 


5 
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station two downstream, Then, any shock would have to be between stations 
one and two. Due to the assumption of uniform properties at station two 
-- necessary for application of the momentum equation between stations 

one and two, where mixing is assumed complete--no definitive information 
can be obtained regarding shock location upstream of station two. It is 
worth noting that should mixing be completed prior to the assumed station 
two, the pressure recovery would be less than that for the situation 
considered in the analysis. For the purposes of this analysis, then, the 
supersonic solution was considered non-existent if possibility (1.) did 
not provide sufficient pressure recovery. 


The well-known normal shock relations, written for a shock at station 





(мо аге: 
P,  2ΥΜΙ-(Υ-ἡ 
τ ΒΕ, 2 (132) 
uie _ 
τ. эь ( + тү м2) (27-м; = ) (133) 
a r 0) 
2 M + = 
M, = z—— LL (134) 
тої M, -I 
Tri ‚= omm 
fe ( 2 M;)T- (rh n) ца (135) 
Pt, T ма +1) т Y dy | 
р, тэг ги МІ — (v nov 
gc s kid. πα i o 
гн 2 vi) | Υ +! ge) 


Here, the subscripts, 72% апа ! ) ', refer to upstream and downstream of 


the shock respectively. For subsonic compression after the shock at two, 


iJ - = 
Аў р, 9 іі 


there ise 


(137) 
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Again, this equation must be solved by iteration to yield p3 and the value 
obtained is compared with atmospheric pressure, As noted above, the solu- 
tion--no shock from station two downstream--is known if ра Ра» If P3 = Pas 
the solution is known and the properties of the fluid at station three are 
determined from (129) - (131), with appropriate subscripts. 

For p3> Pa obtained from the shock-at-station-two-check, the next 
check made is for conditions in the exit plane with supersonic expansion 


to that station, The expansion equation is written as: 
ne t 
и ГЛГЭГ τ == 20 ως Р, ) 4 (138) 
А, Pe, ; ШЕ. А. Pra! 


Here, Pto is a fictitious stagnation pressure determined for stagnation 





conditions which would exist for a converging-diverging nozzle with a 


Mach Number, Mo, at the point where the area is Ay. Then: 
= Пе 
Ne~] 
Pe = P(e)” (139) 


With p} determined from iteration of (138), T4 is obtained from (109), V; 


from (110), written with appropriate subscripts, and stagnation pressure 


from: 
et 
T-i 
і з ВЗ) (140) 
Also, there is: 
Ne „ы (141) 


Гт9КТ | 
The value of P3 obtained corresponds to the value for full expansion, If 
P3 = Pa, the solution is known and the properties and velocity have been 
obtained. ІҒ рз у> Ρα; the flow is underexpanded and the flow adjustment 
required takes place outside the exit plane in the form of a Prandtl-Meyer 


expansion. In this instance also, the solution is known and the properties 
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at the exit are taken as those for full expansion, with comment that 

the noted expansion occurs, For P3< Pa» the practical consideration, two 
possibilities must be considered: an oblique shock at the exit; and a 
shock in the diverging portion. The highest pressure for occurrence of 
an oblique shock at the πρ. is obtained from (132), with the up- 
stream Mach Number equal to that for full expansion. If the pressure 
downstream of the shock is less than or equal to ра, the solution is known 
as a shock at the exit. Exit conditions are then taken as those for full 
expansion, but with comment that an oblique shock occurs outside the exit 
plane. 

If Pg 5 X B4» it is necessary to determine the location of the normal 
shock in the diverging portion. The process of solution is one of itera- 
tion, using supersonic expansion to an assumed distance X from station two, 
followed by a normal shock and subsequent compression. All equations 
required have been given previously, except the one used to determine the 


area at X? ?, From the geometry of the diverging portion there is: 


2 
А, τ η А (r, - r) (142) 


The iteration, using (138), (132) - (136) and (137) in that order, and 
with appropriate subscripts, is continued until the distance X is found 
such that p3 = pas Exit plane conditions are then determined using (129) 
- (131), with the subscript ° 0 * substituted for two. 
Computer Program 

Program EJECTOR, included in Appendix IV, performs the analysis 
described in the foregoing paragraphs in the manner and in the order 
given, Inputs to the program are the values of polytropic efficiencies 
assumed, the number of different inlet conditions for which solution is 


desired, and the sets of values of inlet total temperature (?R) and total 
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pressure (psia). The program obtains the subsonic and supersonic solutions 
for both stator profiles (i.e., the converging and the converging-diverging 
nozzle types) for each set of inlet conditions. Тһе several subroutines 
are listed below in the order called, and the function of each is given, 

l. NOZZLE. This subroutine determines the conditions at the exit 
of the exhauster nozzle by iteration, using (111), (112), (109) and (110). 

2, TURBINE. For an assumed value of pressure at station lp for the 
flow from the turbine hood, this subroutine solves (107) for Pta by itera- 
tion, using also (98), (99) and (100) to determine required flow rates (or, 
flow functions, $ ). Eq. (97) is used to find the total flow rate from 
the turbine hood to the exhauster piping. 

3, MMENTUM. Eq. (125) is solved in this subroutine for the subsonic 
conditions at station two. Im order to get values of velocity and density 
at this station, for use in determining average velocity and density re- 
quired for calculation of the pressure loss, (125) is solved for zero 
friction, Pressure loss is then determined using (115), and (125) is again 
solved for the new value of Pie Eqs, (117), (119), (126) and (127) then 
yield the additional properties at station two. 

4. SUBCOMP. This subroutine solves (128) to determine exit pressure, 
The main program contains the required check for p3 * pa, assumes a new 
value of Ply êS required, and again calls TURBINE, MMENTUM AND SUBCOMP. 
This sequence of events continues until the required condition is met. 

The subsonic conditions at station three are then determined from (129) - 
(131). 

5, ММЕМТОМ, After the subsonic solution is known, this subroutine 
“is again called and the supersonic solution at station two is determined, 
6. SHOCK, This subroutine solves the shock equations for the super- 


sonic solution properties at station two--(132) - (136). 
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7. SUBCOMP. This subroutine is again called for the properties 
determined by SHOCK and the pressure, p3, is determined, If P3<Pa, the 
solution is known as no shock at station two or in the diverging portion, 
If P3 = Ра» the solution is known also and the conditions at station three 
are those given by (129) - (131). In the event that P32 Pa, the program 
proceeds through the subsequent steps, 

8, SUPEXP, Using the MMENTUM solution conditions at station two, 
this subroutine determines conditions at station three for full expansion, 
using (139), (138), (109), (110), (140) and (141), The main program then 
checks Ғог рз = pas If the equality is satisfied, or if p3) pg, the 
solution is known as full expansicn or under-expansion and the exit plane 
conditions are taken as those for full expansion. The output includes a 
Hollerith print of the condition found to identify the type solution 
obtained, i.e., either full:or under-expansion, For pJ < Pa ° the subsequent 
steps are taken, | 

9, SHOCK. To determine the maximum back pressure for a shock at 
the exit plane, this subroutine is again called for the conditions found 
for full expansion to station three, The main program then checks post- 
shock pressure against ambient pressure, The solution is known as an 
oblique shock at the exit if this post-shock pressure is equal to or 
greater than ambient pressure, Such solution is specified by a Hollerith 
statement and the exit plane conditions are taken as those upstream of 
the shock. If the post-sheck pressure is less than pą, the program 
continues with the following calls. 

lO. SUPEXP. This subroutine is called for an assumed shock location, 
X, and with an area given by (142), to determine conditions for supersonic 
expansion to that location, 


ll. SHOCK. Conditions downstream, for a shock at X, are determined, 
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12. SUBCOMP. This subroutine is called next to find the exit plane 
pressure for the conditions downstream of the shock at X, Comparison of 
this pressure with ambient pressure is accomplished in the main program 
as before, If рз = рд, the solution is known, the shock occurs at X, and 
conditions at station three are given by (129) - (131). If P3> Pap the 
distance X is increased and ії рз (Ра» the distance assumed is decreased. 
Iteration varying X continues until the required equality is obtained. 


It must be noted that subroutine NOZZLE is called only once in the 


above series of steps. Hence, it is tacitly assumed that any slight 





expansion or shock at the exit plane of the nozzle is insignificant. 
Additionally, the pressure at station one is assumed constant and equal 
to the pressure for the flow from the turbine hood, since this pressure 
is used in the solution of the momentum equation. 
Results and Discussion 

Solutions of the exhauster analysis were obtained for inlet total 
pressures from 44 through 36 psia in two psia increments, and inlet total 
temperature was varied between 660 and 760°R in 25° steps for each value 
of inlet total pressure. The values of polytropic efficiencies assumed 
are: for the turbine, Ny = 0.75; from the turbine hood to station iw 
Ία-Ξ 0.925; for flow to the nozzle throat, ly = 0.975; from the nozzle io | 


throat to the nozzle exit, Ny = 0.9; for supersonic expansion in the | 


=. 


— «аура 





diverging portion, Ne = 0.8; and for subsonic compression in the diverging 
portion, Ne = 0.7. 
Results of the analysis showed process pressures to be essentially 
independent of inlet total temperature. Of primary interest is the vari- 
ation of total to total pressure ratio across the turbine with inlet total 


pressure, Such variation for both the converging and the converging- 


diverging nozzle stator profiles is shown in Fig. 32. 
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Values of velocities and temperatures given by the analysis varied 
with both inlet total pressure and total temperature, However, values 
of supersonic Mach Numbers obtained for station two, the inlet to the 
diverging portion, were independent of inlet total temperature, Pressure 
in the exit plane for a shock at station two (maximum pressure recovery) 
was also independent of inlet total temperature. Maximum values of P3 of 
13.88 and 13.83 psia for the converging and converging-diverging guide 
vanes, respectively, were obtained for the maximum inlet total pressure 
of 44 psia . Hence, the analysis performed indicates that any shock 
adjustment will occur in the constant area duct upstream of station two. 
This result would seem to substantiate the expectations from general two- 
dimensional considerations that a series of oblique shocks followed by a 
normal shock will occur in the duct, 

Values of referred flow rates, w/6/6 ; obtained from the analysis 
were: 3,535 lbm/sec through the nozzle; 1.438 lbm/sec through the turbine 
for the converging guide vanes; and 1.022 lbm/sec through the turbine for 
the converging-diverging guide уапе5, The totals for the two types of 
guide vanes, in the order given above, are 4,973 and 4,557 lbm/sec. 

Although at present an upper limit of 6609R exists for turbine inlet 
total temperature, during future tests with the exhauster operating it 
may be found desirable to increase that limit to prevent undue condensation 
and/or freezing in the turbine discharge flow. Consequently, provided in 
Fig. 33 are curves of minimum inlet total temperatures required to main- 
tain turbine hood temperature ten degrees above freezing and to maintain 
turbine hood flow temperature at the nozzle exit plane at 492°?R, These 
values are plotted versus inlet total pressure, The curve for Тіт = 32°F 


is common to both sets of inlet guide vanes, 


$5 


5. Determination of Flow Rates. | 
General 

In Section 2, the several required flow rates were enumerated. In 
addition, the method for determining the flow rate through the exhauster 
nozzle was specified (see p. 21). To estimate this flow rate, (111) given 
in the preceding section can be used. In this section, the equations for 
use in determining the flow rates through the flow measurement nozzle, the 
plenum labyrinth seals and the shaft labyrinth seals are developed,  Equa- 
tions developed for the flow nozzle and plenum labyrinth leak rates are 
based on experiments, but the equation for shaft leakage provides flow 
rate estimation only. 

Test Procedures 

The configuration for flow nozzle calibration tests is indicated in 
Fig. l, p. 21. The removable standard orifice assembly was installed, 
the outlet from tank number two to the turbine and the two-inch pipe were 
blocked, and the valve for controlling the flow to the exhauster nozzle 
was closed, Тһе standard square-edged orifice of the removable assembly 
is fitted with both vena-contracta and flange taps, Upstream pressures 
were measured against atmospheric reference in inches of mercury, while 
differential pressures were measured in inches of water using a portable 
manometer with a range of 70 in.,and with 0.1 in. graduations. Flow 
nozzle instrumentation is shown in Fig. 3, p. 24. The accuracy of 
pressures and differential pressures was taken as + 0.03 in. and the 
accuracy of temperature measurements was assumed to be +1.0°R. 

Flow nozzle calibration tests were made at essentially constant 
supply pressures of either 30 or 45 psia. Flow rate was controlled by 
using the valve at the exit of the removable standard orifice piping to 


set the differential pressure accross the vena-contracta taps. Supply 
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pressure was maintained by adjusting the main discharge valve of the 
compressor installation. Differential ргеззигей was increased in 
increments of ten inches of water for increasing flow rates to the 
maximum differential pressure obtainable, then decreased by five inches 
of water, followed by incremental reductions of ten inches of water to 
minimum flow. Approximate ranges of flow rates of the tests were 1.44 to 
3.67 lbm/sec at 30 psia, and 1.8 to 4.6 lbm/sec at 45 psia. 

Instrumentation of the two-inch pipe orifice, which is used for deter- 
mining plenum labyrinth seal leak rates, is basically the same as that of 
the removable standard orifice assembly. Plenum labyrinth leak rate 
tests were made by varying the supply pressure (approximately) from 24 to 
44 psia. Pressures were set by observing the turbine plenum reading. 
Supply pressure was varied in increments of about five inches of mercury 
for increasing pressure; then, after an initial reduction of 2.5 in, Hg 
from the maximum supply pressure, incremental decreases of five inches of 
mercury were used, Tests were made without the exhauster and associated 
piping installed. Hence, the maximum ratio of seal inlet total to static 
discharge pressure was limited to approximately 3. 

Formula Development 
The basic flow equation for a square-edged orifice or flow nozzle, 


in lbm/hr, 151 


= 2 
W- 359.1 D. «КҮ, ^ hw (143) 
Here; 
D, = orifice diameter or nozzle throat diameter (in). 
œ = coefficient of thermal expansion based on р2 (dimensionless). 


"Stearnes, R. F., et. al., Flow Measurement with Orifice Meters 
(New York: Van Nostrand Co., Inc., 1951) p. 6. 
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discharge coefficient (dimensionless). 


ж 
Il 


Yq] = expansion coefficient (dimensionless). 
PS fluid density at upstream tap (1bm/ft?). 
hy = differential pressure across taps (in. H,0). 


The equation of state for an ideal gas is: 


e 0) (144) 


Неге, ру is in psia, R - 53.35 ft-lb/lbm-^R, and T, is in °R., For pj in 


inches of mercury, (144) becomes: 


о Р (о.4912)(144) _ й, 
79^ ΠΠ τα δες (are) 


Inserting (144a) in (143) and converting to lbm/sec: 





W= 0.115 D xK Y, Ым, (145) 


The coefficient, Cx, is determined from the linear coefficient of thermal 


? e ; 
expansion, О,. For р, = orifice diameter at elevated temperature: 


13 


Neglecting the higher order final term and letting (22)? т ос 02 : 


£ 2 
= [I+ AÐ] =| + 20% (At) + (At) ay (146) 


x D; =[ | +2% (4t)] 0 (147) 
From (147) it is evident that; 

X = | 2 % (At) (148) 
For both standard orifices, made of type 304 stainless steel, 


СХ « (9.5)1079 9Е71, giving: 





OÇ = | + 0.0019(5,588) (149) 


(The value, 0.00193 was used ir TN В For the flow nozzle, made 


"Ibid,, p. 257. 
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of 2024-T4 aluminum, © = (12,6) 1079 ?r-l, ana: 
Oy, = 1+ 0.00252 (+588 | (150) 


The discharge coefficient, K, is dependent on the ratio of orifice 
or nozzle diameter to upstream pipe diameter (і.е., 6 - D5/14), Reynolds 
Number based on 02, and the type measuring taps used. To simplify 


computations for square-edged orifices, the following relation may be 





used:? 
K IO 

T+ A/R) ~ U HARE (151) 
where ;1O 

A = D,(830 - 50008 + 900062 — 42090” +B) (152) 
with: 

B - PX for flange taps (153a) 

p == - 100 for vena-cont, taps (153b) 
The Reynolds Number expression can be written: 11 

Re = 630. | (154) 


where zt is absolute viscosity in centipoises at flow temperature and 


pressure, For zt independent of pressure and assuming linear variation 


of zt in the temperature range from 50 to 300°F, there 18:12 
а - |008'- 1.9 - 0.24 ЖЭ En (155) 


Putting (155) in (154), with W - 3600w, and defining, X є В. (1079): 


D, z 


9Ibid., p. 65.  lÜIbid,, p. 212.  irbid., p. 64. — l^Fbid,, p. 335. 
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Let §=1+A/R,. Then: 
-b 
ΠΠ, або (157) 


Eq. (151) then becomes: 





4, 
K K ( 
— = 158) 
5 f^ 
For flange taps, values of K* - Коф corresponding to Re* = 00, are 
tabulated. 13 Eq. (158) for flange taps is then: 
= f Ko (158a) 


For vena-contracta taps, this shortcut method is not authorized and values 
ої Кор are not tabulated., With K tabulated as a function of Reynolds 


Number, 14 


an iterative procedure using these variables is normally 
required, with flow rate, w, as a third variable. However, by judicious 


choice of Re”, and hence K*, (148) becomes for vena-contracta taps: 


E. uS (158b) 
For the removable standard orifice assembly, Re* was chosen as (105), 
giving K* = 0,698. For a Reynolds Number of 2(10°), the value of K deter- 
mined from (158b) differs by 0.114 percent from the tabulated value, while 
the range of Reynolds Numbers for tests was approximately 3.9(10?) to 
1.3(106). Similarly for the two-inch pipe square-edged orifice, Re* was 
selected4(105), giving K* = 0.6103. At a Reynolds Number of 5(104), the 
computed value of K differs from the tabulated value by about 0.1 percent, 
and at Re = 2(10°) the difference is about 0.13 percent. The range of 
Reynolds Numbers in the tests conducted to date was approximately 5(104) 


το 1.2(105), Tests yet to be conducted with the exhauster installed will 


14 44,, там 20. 


l3Ibid., Table 19, 
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be at somewhat higher Reynolds Numbers than 1.2(10?), However, since it 
is anticipated that the labyrinth seal flow will become choked at pressure 
ratios (plenum total to seal discharge static) only slightly higher than 
3:1, Ве з 2(105) is considered a realistic upper limit. 

The expansion factor, Үү, compensates for compressibility in the 


flow through the orifice or nozzle and is defined as;15 


Y, = 1.0 — (0.41 + 0.358") =P (159) 


Here, Ap and p, must be in consistent units, and Y is the specific heat 
ratio, For T= 1.4, Ap = h, (in. H,0), and р) з Pi (in. Hg), this 


equation becomes: 


4 
(-10-44 00:80 (160) 


(1.4) (13.59) 


Eq. (145) can be written in the form: 


М = Co $Y |Қ» (161) 


Here, the constant, C, includes р2 and either Ко , for flange taps, or 


E. for vena contracta taps. Values of C, along with values for A of 
(152), and the several constants in the expressions for Yj and X, were 
determined by using the CDC 1604 computer, Table III lists the values to 
be used for each standard orifice, In solving (161), it is necessary to 
assume initially that 5 - 1. With the resulting value of Ww, X is deter- 
mined from (156), and the required value of Í is obtained from (157). 
Eq. (161) is then solved using this new value of €f. 

Since the object of the flow nozzle calibration was the determination 


of the discharge coefficient, K,, (145) is written: 


& [n 
Kn = 0.5 D? c Yi [P. h. 202) 


l5rbid., p. 71. 
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TABLE III | 
FLOW EQUATION CONSTANTS - SQUARE-EDGED ORIFICES Y 
TRANSONIO TURBINE TEST RIG 


Flange Taps Vena Cont Taps ΠΡΙ 


σι ο τω ° 
ma 
0.026 


EQUATIONS 












A 


б) — $*sox$yx [м | 
Mic 
1 


(152) -A = Da(830 ~ 500984 90004? - 420047 + B) 


| "E : 
s a ü - Зо. P". um 
; = 4 1 | › 
» S | B CUM η h 
160). Ү,:1.0- 0.41 + o. = s «О = ν 
(1 y ice loa (15:57 | Р 1.0 = c; т 


(1555) | 221.94 0.24 Επ 
Me 100 | 


(149) Az 1.0 + 0.00193 ΝῊ 


(157)  : 521.4 _A (107°): 
\ | X 


і 
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Then, with the flow rate, w, determined from (161), οἱ, given by (150), 


and ү, by (159a), K. is found from (162). 


n 
The basic method used for estimating leak rates through labyrinth 

seals is that developed by Egli, [7] The equation of continuity combined 

with the energy equation for an ideal labyrinth seal (i.e., all kinetic 


energy is destroyed after each throttling) is: 





й «оф s... т (163) 
where: 
: 40 2 % 

Фе ELA. |= aac] эн" 
aX = discharge coefficient for a single throttling (dimensionless). 
A = cross sectional area of seal passage (іһ2), 
Реог total pressure at seal entrance (psia). 
p = static pressure at seal discharge (psia). 


Те = total temperature at seal inlet (^R); Tt, = constant across the 
seal for assumed adiabatic flow. 


R = gas constant for air (ft-lb/lbm-?R) 
n - number of throttlings. 
r = overall seal pressure ratio (р/р). 


To account for a non-ideal labyrinth seal, a velocity carryover factor, 
T is introduced to adjust for carryover of kinetic energy from one 


throttling to the next. Eq. (163) then becomes: 


АР, 14 
СЕЕ 6 иа 165 
οί 0 πα R ( ) 


For §= tooth clearance, τ = tooth width, and s = chamber width, (see 
ον απο 25): 
ot = Ғ(6/с, п) = constant 


` 9 ti j = 
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T's f( Š /s, n) = constant 

The dimensionless seal pressure ratio function, (0 , was developed 
by series expansion with a lower limit of pressure ratio across any one 
throttling of O.B, Since preliminary analysis of exhauster performance 
showed an overall pressure ratio across the shaft labyrinth seal of about 
0.55, this assumption was considered valid for the ten throttlings of 
that seal, However, since plenum labyrinth seal leak rates were deter- 
mined experimentally, use of (165) was limited in this instance to 
attempts to correlate experimental data with the method used and to deter- 
mine whether test data for overall seal pressure ratios of 0.34 to 0,6 
could be extrapolated to lower pressure ratios of 0.16 to 0.34. For the 
two identical labyrinths on either side of the plenum, (165) can be 


written, introducing the non-dimensional flow function, ? : 


м Л [К 
ф = Arp, 4 = К, (166) 


If x and т“ аге іп fact constant, K - constant, and experimental data 
for higher values of r can be extrapolated to lower values of r. Since 
it is anticipated that the flow through the seal will be choked at 
( яах = 0.2798 (r = 0.225), the range of extrapolation required would be 
relatively small. 

The curves of c versus $/t given by Egli [7] make determination of 
C for the geometry of the shaft seal somewhat uncertain, Consequently, 
the value of @X = 0,76 was taken from experimental data presented by 
Jerie, [8] where & is plotted as a function of $/t and 6/1, Неге, 
b = chamber heighth (see Fig. 4, p. 25). Plots of T* versus ó/s given 
by Egli are based on experimental data for $50.01 in, Egli'!s citation 


of Friedrich’s tests of straight-through labyrinths with = 0.006 to 
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0.01 in,, and with two and three throttlings, indicates that somewhat 
smaller values of  Y* than those given by Egli would be proper for the 
geometry of the shaft seal of the TTTR. Considering also, ten rather 
than two or three throttlings, and data for multiple throttlings given 

by Jerie, Friedrich's value of 1.15 was selected for Y* in preference to 
Egli's value of 1.25. With these values of c and Y*, the final formula 
for estimating shaft seal leakage is, with A = T (1.2)(0.005): 


% = 00128 Ав (9, (167) 


| Teo 


Then, introducing the non-dimensional flow function, 9 : 


9 T 7 
φ. - жыла = 0.874 Q; ; T, E ТЕ 5 а A (168) 
to 


Since Y= f(r), the shaft seal leakage can be represented by a single 
plot of $ versus overall seal pressure ratio, r, This plot is show in 
Fig. 34. 

Results 

Since the tolerance for vena-contracta taps is less than that for 
flange taps, plotted and tabulated data included herein are vena-contracta 
values unless otherwise indicated. 

Fig. 35 shows nozzle discharge coefficient, K,, versus Reynolds 
Number as determined from tests of 8 September and 13 September 1965, It 
is evident from this plot that K, can be taken as constant at Reynolds 
Numbers greater than about 6(105)--(% € 2 lbm/sec). Values of K, = 
constant were determined by the method of least squares for the following 
data groups: 8 September tests at 30 psia; 13 September tests at 30 psia; 
these two tests combined; 13 September tests at 45 psia; and all listed 
tests combined. Resulting minimum least-square-error values of К, аге 


tabulated in Table IV. Im addition, values of Kp = constant for vena- 
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Pressure Ratio,r 


SHAFT SEAL FLOW FUNCTION 


FIGURE 34 
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Reynolds Number, Re(106) 
FLOW NOZZLE DISCHARGE COEFFICIENT 


TRANSONIC TURBINE TEST RIG 


FIGURE 35 





TABLE IV 





Flow Nozzle Discharge Coefficients 


Transonic Turbine Test Rig 


Date Pressure | Points| Error #| | | Ка | Error £ K 











| | .30 psia | 12 | .01055 | 1.0026 | ,00684 | 1.0012* 
30 psia 21 201531 1.0012 .01071 1,0024 
30u pai a 21 __ |. — a ee ee 
8/13 Sep 45 psia 12 ‚01614 1,0011 ‚01633 1.0012 
13 Sep 45 psia 12 ,00497 1.0008 .00864 2 а9992% 


oe 


# Error figures show least squares average difference of data points 
from Kn for minimum error 


* Recommended values 
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contracta data for the individual tests are indicated in Fig. 35. Recom- 
mended values of K, for 30 psia and 45 psia are indicated by an asterisk 

in Table IV. Some difficulty was encountered with water collecting in 
orifice pressure tap lines during the 8 September tests (30 psia). Although 
these lines were continually drained during such tests, the data are 
considered inferior to those obtained at 30 psia on 13 September, when 

the water collection problem did not exist. 

Experimentally determined values of K, and ф, for the plenum 
labyrinth seal (see (166)) are plotted versus overall seal pressure ratio 
in Figs. 36 and 37 respectively. It is apparent from Fig. 36 that extra- 
polation of K, to lower pressure ratios would be of questionable accuracy. 
Hence, subsequent tests with the exhauster operating will be required. 

To provide a mathematical expression for Ó, = f (r), the data plotted in 
Fig, 37 were fit with a parabola by the method of least squares, The 


general equation used in the curve fitting was: 


і 


In addition to experimental data, the known condition that at r = 1,0, 


- 6546 (169) 


ф, з 0, маз introduced. Since all values of $, were of like sim (+), 
symmetry about the r axis was assumed by setting b = O, The tests of 
20 October 1965 showed considerable variation between values of K, and $ 
computed from data for increasing supply pressure and decreasing supply 
pressure, Consequently, ¿special care was taken to assure flow stabili- 
zation prior to recording data during the 26 October tests. Although 
the noted discrepancy was not eliminated, it was reduced materially. 
Hence, the 20 October, decreasing supply pressure values of 9, (corres- 
ponding to K, Points to the upper right in Fig. 36) were neglected in 


the curve fitting process. The resulting least-squares second order 
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Pressure Ratio, r 


LEAST SQUARES PARABOLA 


Seal 


FOR THE 


PLENUM LABRYINTH SEAL FLOW FUNCTION 


FIGURE 37 
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polynomial obtained, for which the least squares error is 0.02079, is: 


2 


Solved for the normally dependent variable, Ó, » (170) becomes; 


αι ο στ 
А 4.0191 (170а) 


This equation is plotted in Fig. 37 along with the experimental data. 
Also shown in Fig. 37 are segments of three curves, ᾧ,- 2κ1ρ, for values 
ої К, appropriate to three different ranges of seal pressure ratio, 

Data reduction required and the solution of (164) were accomplished 
by using the CDC 1604 computer, Listed below with their purposes described 
are the various programs used, Тһе programs themselves are included in 
Appendix V, along with the description of inputs and outputs for each 
program, Since these programs use the equations of this section in a 
straight-forward manner, no detailed explanation is given, 

1, ORIFICE. Computes the constant C of (161), the value of A of 
(152) and the constants in the expressions for Y, and X. 

2. FLOCAL. Determines the values of K, of (162), based on experi- 
mental data. 

3. FLOCAL2. Finds the least squares value of K, = constant, 

4. LABRINT. Solves (164) for pressure ratios from 0,025 to 0.975. 

5, LABLEAK. Determines flow rates and values of Ky, and $, of 
(166), using experimental data. 

6. MINSQ2. Determines the least squares parabola for plenum 
labyrinth leakage (169). 
Discussion and Recommendations 

Based on the results of nozzle flow rate experiments conducted, it 

is recommended that a value of 0.9992 be used for the nozzle discharge 


coefficient in (145). However, for supply pressures of 30 psia or lower, 


111. 


a value of 1.0012 would be more appropriate, 

Eq. (170a) is recommended for determination of plenum labyrinth 
leak rates in the pressure ratio range: r = 0.03 to 1.0. Since plenum 
labyrinth leak rates (0.04 to 0.1] lbm/sec) are very much smaller than 
turbine flow rates (2.0 to 3.9 lbm/sec), the slight errors in leakage 
flow rates introduced by the use of this equation should have negligible 
effect on the accuracy of turbine flow rates. Furthermore, this equation 
should give satisfactory flow rates at pressure ratios, r, higher than 
those used in its generation. 

For exhauster operation with a supply pressure of 45 psia to the 
exhauster nozzle and the turbine, the estimated shaft seal leak rate 
given by (167) is about 0.003 lbm/sec, while flow rates through the 
exhauster nozzle and the turbine are estimated as 9,38 and 3.87 lbm/sec, 
respectively. Using Eglits [7] values of discharge coeffioight and 
carryover factor ( “= 0.9 and “ғ z 1,25), shaft seal leak rate estimates 
would be about 50 percent greater than estimates given by (167). However, 
plenum labyrinth leak rates predicted using values of X and γ᾽ deter- 
mined from Egli are approximately 50 to 70 percent greater than the leak 
rates determine experimentally., Hence, (167) is recommended for predicting 
shaft seal leakage flow rates. 

The recommended flow rate equations are given in Table V for ready 
reference, 

Tolerance in Flow Measurement 

The overall error tolerance in percent in flow measurements using 
бай edged orifices is given by the square root of the sum of the squares 
of "exponents times tolerances, in percent" ascribable to the individual 


items in the flow equation. Непсе:16 


16ть44,, р. 24. 
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TABLE V 


SUMMARY OF FORMULAS FOR DETERMINING FLOW RATES 
TRANSONIC TURBINE TEST RIG 


| 
Flow Rate. Basic Formula : Auxiliary Formulas 
Nozzle ! і F τ 
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д = (L24 (0 + [A(e0] “(АК ЕГАР 


АЫ Ην (171) 
Tolerances of the individual items for the Transonic Turbine Test Rig 
installation are given in Table VI, 

The error in values of discharge coefficients of the flow nozzle and 
the square-edged orifice must be approximately the same, It follows that 
the tolerances in flow nozzle flow rate measurements must be the same as 
those for the square-edged orifice used in calibration. Hence, vena- 
contracta-based values are preferable, due to a lower error tolerance 
than that for flange taps. 

Tolerances in measured labyrinth leakage rates also should be de- 
pendent on the tolerance of the calibration orifice. However, it is 
apparent from Fig. 35 that there is considerable scatter in the data, 
Hence a tolerance arbitrarily set at + 10% would appear reasonable, 

Error tolerances in estimated shaft seal leak rates are difficult 
to specify exactly. Here, a reasonable tolerance could be set at + 20%, 
based on the qualifications set forth in the literature cited, 

6, Calibration of the Force and Moment Measurement System 

The axial force and moment about the turbine axis of the stator- 
plenum assembly are measured by Wiancko, two-arm, variable reluctance 
bridge force pickups using a phase shift frequency modulated oscillator. 
The linear range of the oscillator is from ten to 12kc. The permissible 
loadings of the two capsules are 50 lb, and 600 in-lb for the force and 
moment pickups, respectively. 

Tests were conducted to establish the linearity conformance of the 
measurement systems and to determine whether or not cross-coupling exists 


between the two systems. In addition, readings of force and moment were 
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TABLE VI 


PERCENT TOLERANCES OF FLOW EQUATION VARIABLES 
AND PARAMETERS - SQUARE-EDGED ORIFICES 
| TRANSONIC TURBINE TEST RIG 
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Independent of the Orifices . 


A(X) = £ 0.05 : s 
| hy 


А(1,) =-0 tot0.05 Рог 0.0<— < 0.015 
| Р, 
з 20,5 toil.O | for 0.015 « b. 
. h Pi 
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A(h,) = 9:92... (100) 


ο £8 
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Removable Standard Orifice Assembly: 
= + 0.0005 \ 
A(D2) = + 52425 (100) = + 0.0118 
A(X) EI L0 for flange taps 
|  — 
z £ 0.4 to = 0.6 for vena contracta taps 
Two-inch Pipe Orifice: 2 


A(D2) = + 0.0002 (100) = + 0.0606 
|". 085 t | 


.. AK) -». X 1.0 | for flange taps 


\ | 
з ® 0.4 to X O.6 | for vena contracta taps 
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taken during plenum labyrinth leak tests to check for net force and/or 
moment due to the flow through the plenum labyrinth seals, 

Fig. 38 shows schematically the configurations for the force and 
moment calibration tests. Standard weights were used to obtain known 
forces and moments, The following description of procedures applies to 
the force system; however, the same procedure is used for moment calibra- 
tion, except that maximum load is 600 in-lb and increments of load are 
100 in-lb rather than the corresponding values of 50 lb and 12.5 lb for 
force capsule calibration, The zero setting was obtained with no load 
as ten kc and, with the full load of 50 lb, the bandwidth was adjusted 
to set the reading at l2kc. Linearity was checked with one-half full 
load applied, with the desired reading as llkc. If this reading was not 
obtained, the linearity was adjusted so that the new reading was differ- 
ent from llkc by three times the original difference from llkc, and in 
the same direction as that difference. The zero point, bandwidth and 
linearity were then adjusted in order repeatedly until the desired read- 
ings were obtained without further adjustment (i.e., ten, 12 and llkc 
for zero, full and one-half loads, respectively). Starting with zero 
load, the force was increased in increments of 125% со the maximum load 
and then decreased by the same increments to zero load, with readings 
rotos for each load. The differences in readings for loading and 
unloading at each given force indicate the hysteresis (and non-linearity) 
of the system, The maximum departure observed for the force capsule was 
eight cycles per second (0.2 1b) at the one-half load point; maximum 
departure for the moment system was five cycles per second (1.5 in-lb). 

Cross-coupling between the two systems was checked by applying a 
fixed load to one system and thenapplying various loads to the other, 


while watching for changes in reading of the former, Although it was 
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concluded that no significant cross-coupling exists, there remained some 
doubt as to whether the moment affected the force reading. The uncertain- 
ty arose because unpredictable changes in force reading were noted, These 
changes in reading were consistent neither in magnitude nor direction; 
also, by jiggling the moment calibration arm the reading could often be 
returned to the proper value, It was concluded that the variation of 
force readings was caused by the indeterminate friction force developed 
between the shaft to which the moment calibration arm is attached and the 
bearings through which the shaft passes, It is significant in this regard 
that full-load displacement of the capsules used is only 0.05 in, 
Pertinent also is the fact that no changes could be detected in the force 
reading while applying moments with a wrench, and with the aforementioned 
shaft disengaged. 

During plenum labyrinth leak test, no significant changes in either 
force or moment readings were observed. It was concluded, therefore, 
that the flow through these seals produces no net force or moment. 

6, Turbine Tests. 
General 

Tests of the turbine were conducted both with and without the rotor 
installed. Only the converging nozzle profile stator has been tested and 
tests with the rotor installed were done with this same stator and the 
full-size, uncut rotor, Axial clearance for tests conducted to date was 
0.243 in. The gas path for these tests was the basic path shown in Fig. 
10, p. 34. Rotor tip clearance was 0.025 in. Instrumentation for the 
tests was as shown in Fig. 6, p. 28, except that the total pressure probe 
line apparently became disengaged inside the machine, Hence, subsequent 
to the third rotational test, the total pressure upstream of the stator, 


Pty» was not measured. Values of Pt, were determined by applying a 
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correction, based on previous runs, to the static pressure at the probé 
location, pg. It bears repeating that the pressure on the downstream 
side of the cover plate inside the stator was measured during stator tests 
but not during tests with the rotor installed. 

Prior to each test, the force and moment measurement systems were 
calibrated for proper readings at zero, full and one-half loads, as 
described in the preceding section, Also, prior to each test with the 
rotor installed, the rotor moment measuring system was calibrated for: 
minimum reading of zero at 50 in-1bs; and maximum reading 100 at 700 іп-1ь, 
Then, the applied moment was increased from 50 to 700 and then decreased 
to 50 (in-lb). Readings obtained were plotted versus applied moment to 
yield the calibration curve. During all tests it was necessary to main— 
tain a tare loading in the axial direction to keep force readings in the 
linear range, 

During stator tests, inlet total pressure was varied from minimum 
values to the maximum achievable ahead of the stator, about 42.5 psia. 
Data were recorded at increments of supply pressure of about five inches 
of mercury. 

Since it was planned to present turbine performance data in the form 
of curves for constant pressure ratio (total inlet to static discharge), 
for tests with the rotor installed the dynamometer was set initially for 
maximum power absorption, Supply pressure was increased until the value 
for the desired test pressure ratio was obtained. With desired pressure 
set, and with maximum power absorption, minimum RPM for the test pressure 
ratio resulted, Then, RPM was increased by reducing power absorption and 
data was obtained at desired values of RPM, while turbine inlet total 
pressure was maintained essentially constant, Inlet total temperatures 


were kept below 200^F, while the minimum inlet total temperature used-- 
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computed in advance--was that value for a static temperature of 42°F at 
the rotor discharge. 

Experience showed that supply air pressure and temperature and RPM 
each affected the others to a considerable extent, Hence, every effort 
was made to record RPM, rotor moment, stator force and moment, and the 
several pressures and temperatures as simultaneously as possible. 

Data Reduction 

The turbine flow rate is determined by subtracting the plenum laby- 
rinth leak rate from the flow rate through the flow measurement nozzle. 
The former is obtained from (170a) and the latter from (145), with K,, 
the nozzle discharge coefficient, equal to 0,9992, 

As stated in Section 2, the instrumentation of the stator-plenum 
assembly was designed to permit use of the momentum and the moment of 
momentum equations to determine stator performance. In vector form, the 


general momentum equation, for the assumed steady flow, is: 


KL -| dm. = [-5 gas, T [Fen ds, + (172) 
Similarly, the moment of momentum equation can be written: 


[^ x V, dm, — f, x Vp dine = [хт p, ds, +] x Tie P, dSe + M (173) 


The integral term for shear forces between fluid particles in the exit 
plane of the stator has been omitted from each of these equations., Such 
forces act perpendicular to the turbine axis and would have no effect in 
the momentum equation for the axial direction, Also, their effects in the 
moment of momentum equation would be small, since the forces at either 
side of the between-blade flow channels oppose each other and should be 
nearly equal, In (172) and (173), the terms F and M represent the force 


and moment, respectively, exerted on the flow by the blades and walls, and 
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are due to pressure and shear forces, 

The control volume used consists of the inner flow passages for the 
fluid, with openings at the inlet, station te’, and at the stator exit, 
station one, Fig. 39 shows these stations, along with the several pres- 
sures, dimensions and forces pertinent to the reducing (172) and (173) 
to the forms used in data reduction, Fig. 6, p.28, shows the stator- 
plenum assembly in more detail. 

It is evident from Fig. 39 that integrals for the entrance station 
contribute no axial force, since the pressure acts perpendicular to the 
axial direction, and since the velocity of the fluid entering the control 
volume has no axial component; too, provided that the pressure and velocity 
at the inlet opening are symmetrical with respect to the turbine center- 
line, no net moment about the centerline acts on the flow at the inlet. 
Hence, the integrals containing pressure and velocity at the inlet 
reduce to zero in the equation for the moment of interest also. The 


expression for the axial force then becomes: 
F= Va, Ч Р, A, (174) 
and that for the moment about the centerline reduces to: 


Vu 


M,- 15 (175) 


Here, V 


q, and V, are the values at the mean radius, гу. 
l m 


l 
The force and moment pickups, affixed as shown in Fig. 39, measure 
the reaction force and moment acting on the stator-plenum assembly, 


Hence, these pickups measure the negative of F, and Μα. Letting R; and 


M; denote the reaction force and moment, (174) and (175) become; 


ο. UD (176) 
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Μι ET mn (177) 


Included in the force measured by the axial force capsule is any 
differential force acting on the exterior of the assembly. As shown in 
Fig. 39, ambient pressure acts to the right over the entire case, while 
acting in the negative direction are the pressure on the cover plate, 
pir? the pressure pj,, acting on the area (Ac-Ap); and, presumably, 
ambient pressure acting on the annular area (A,’ -A,). Let R be the 
reaction force measured by the force pickup, taken as positive in the 
positive i, direction. For R positive, the force reading given by the 
force measurement system will be positive (i.e., greater than сеп/Кс). 


Then, since the sum of all forces in the axial direction is zero: 


R* mV, - p Ås + p Åe- Үс А,)- p, A, -0 (178) 
Solved for the desired quantity, V4,, and substituting w/g for in, there 


Va = W LRE pA -BA R A A) p Aj] (179) 


Considering that a positive moment acting on the stator-plenum 
assembly results in a negative reading from the moment measurement system 


(less than ten kc), for M, defined as the measured moment for the stator: 


ED. (180) 
Solved for Vu» and with w/g for m: 


М = - М, (181) 


W Tm 
With an imposed axial load (tare), R in (181) must be the force given by 
the reading minus the tare force, 


With V4, and Vu kuown from the solutions of (179) and (181), the 


1 
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momentum mean flow discharge angle at station one is obtained from: 





c,-tan' y (182) 
Then, there are: 
VS m (183) 
! cos %, 

qs 

ее a 247; (184) 
Y-! 

р. BE 

Ls Т, | 5 і ал 


and, stator efficiency is given by: 


Τ 
ο 
ies cu. (186) 


For choked flow conditions through the stator, an area restriction 
factor, Ё , can be obtained using the non-dimensional flow function, $. 
This rea restriction factor indicates the fraction of the throat area 
which would be required to pass the actual flow rate under isentropic 
conditions and is analogous to the factor, Ё , used in the turbine 
performance analysis (Section 3). This restriction factor is defined as: 


зі | R 
з Knis ἃ 


— = = = 187 
Ë w* Ma rs Ё ῥ᾽ ( ) 
AR, (89 





Хо : : 
Here, ᾧ , is the choked-flow isentropic value of ф for the throat area, 


A and w is the isentropic flow rate for choked flow. Then; 
+ 2 E 2Υ 
Ó = (+47) {Fi (188) 


The non-availability of the throat pressure is the reason why the re- 





striction factor can not be found for non-choked flow, since, for such 


condition, $; s must be found from; 
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In addition to the foregoing, an effective exit plane area was deter- 
mined using the axial velocity given by (179) and continuity considerations. 
The continuity equation, solved for this area is: 


МКТ, 


P V. (190) 


Ae = 


When compared with the geometric area of the exit plane annulus, this area 
gives an approximation for the blockage of flow area due to both boundary 
layer effects and trailing edge thickness. 

To be calculated for turbine tests ares the momentum mean values of 
absolute and relative velocities, angles and properties after the stator 
and rotor; turbine horsepower; and the efficiencies of the stator, rotor 
and the complete turbine, 

Since the pressure, Pli is not measured during tests with the rotor 
installed, use of the momentum equation (179) must be considered approx- 
imate only. For comparison purposes, however, Ply is estimated so that 
Ya) can be obtained from(179) for use in (182). Due to viscosity, the 
fluid between the rotor and the cover plate will have some mean value of 
angular velocity, which can be expressed as a fraction, f, of the rotor 
angular velocity. Any value of f assumed is, of course, quite conjectural, 
However, it is known that the fluid adjacent to the rotor has an angular 
velocity equal to that of the wheel, while the fluid adjacent to the cover 
plate is at rest. Fig. 40 shows a fluid particle of unit depth into the 
paper. The forces acting on the particle are the centrifugal force due 
to the particles angular velocity, Uf= fury , and the pressure force. 


The pressure forces arves 
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PARTICLE RADIAL EQUILIBRIUM 


FIGURE 40 
Е, (Р 3st *E)do - (ps 3E dc - 4f) d6 + 2 p48 ar (191) 


2 


Since the centrifugal force per unit mass is оу 


r, the centrifugal 


force acting on the particle is: 


F 


e 


‚= ө r(erd8dr) (192) 


For radial equilibrium of the particle, the pressure forces must equal the 
centrifugal force. Expanding (191), equating (191) and (192) and dividing 


by rd@ dr: 
de- = eur (193) 


Integrating from r = SE where p * Ply? to r, where the pressure is p: 


Р-8, 7225-08-14) (194) 
ог: 
р авг 7258 (ru — r^) (195) 
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The differential axial force acting on the cover plate is? 
4 
dF, = e dA - [n, - £2. (rà -r*) | zer (196) 


This equation, integrated from r = O to r = rly? and divided by the area 
of the cover plate yields the expression for determining the average 


effective pressure to be used as Ply? 

















то "n = 
Integrated, (197) is: 
Ши ος MID PS ПЫ 
ше e NL ν.μ”, αὐ : 
pues nz (198) 
й 
ог: 
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For Рін in psia, /2 іп lb-sec^/ft^, and m in inches, the factor, 1442 
in, converts the second term of (199) to psia. Expressing Wp as a 


fraction of the rotor angular velocity, (199) becomes: 


=P. (te) r, (200) 


б — ^2 44) 


I 
Using assumed values of f = 0,5 and bis 0.002378, the pressure 
given by (200) is used in (179) to obtain an approximate value of the 
flow angle at the stator exit and the fluid properties at station one, 
However, the values obtained are not used in subsequent calculations. 
Rather, the conditions at one are obtained from continuity considerations 
and the value of Vu given by (181). The isentropic expansion temperature 
at station one is determined from (186). An iteration process is then 


used. The velocity at one is given bys 


У, 2 29ο (l — 1) (201) 
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Values of Τη are assumed by adding increments to Т ів" The velocity is 


determined from (201), and the flow angle at one is given by: 


201011 13 (202) 


Then, a value of πι is computed for comparison with the value assumed: 
1251 


bas . - (204) 


' 
When a match for Tj - T, is obtained, the velocity and temperature are 


known, and V,, is found from: 


a] 
Va, = V, cos &, (205) 


The solution for the relative flow velocities and angles ahead of 


the rotor is obtained using the velocity triangle relationships: 


W. = Мечі, ο “ως μ΄ (206) 

A = tas (ta (207) 
_ _ X. 

W, = Te (208) 


Equivalent enthalpy is determined from: 


2 2,4 
ιο, E Ч "AT = ш, /12 (209) 


With p, known as ambient pressure, conditions after the rotor are 
obtained by the iterative process;-based on continuityg-used for the 


stator, but with appropriate relative flow variables and equivalent 





enthalpy. The required relative peripheral velocity, Wuz? is determined 


from the measured rotor moment, Mp. The change in total enthalpy across | 
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the rotor is given by: 


AH=H,-H, = 629 | (210) 


w 


Here, AH is in ft*/sec*, Then, from Euler?s turbine equation, written 


for positive AH: 


νως Мо 4 (211) 


The required value of Wu is then found from: 


Wu, = Vu, ~ Up (212) 


The isentropic temperature used as a base for assuming values of T, is 
that value, Το.» for isentropic expansion from Тер апа Peg to the pressure 


po. There 15: 
τα] 
Low Нь (213) 


and: 
ΎὙ-ι 


" р\ т 
Τις he, НЫ (214) 
When the assumed and computed values of T, agree within tolerance, 
the relative velocity, flow discharge angle, and temperature at station 


two are known, The following equations are used to complete the solution: 


20007 E (215) 
275 22 (216) 
шаг 


Referred values are obtained, as is done in the performance analysis, 


using: 
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P. 
0 = 725 -—- 
M 
Ше = су (220) 
Мы = re (221) 
ар” 9 (222) 
MEN. 
ΠΝ les (223) 
kis = -Ahi : Аһьіп ес (224) 
апа: 
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The quantity, Ki,» is defined as the head coefficient and is related to 


the often-used velocity ratio, U,/Cy; where c2/2 = ЛЫ. Бу: 


| | 
t ышар | (226) 





The theoretical degree of reaction, r*, is the ratio of isentropic enthalpy 
drop from p to po (determined along the line of constant entropy given by 
stator inlet conditions) to the isentropic enthalpy drop from Ре to P2e 

ο 


The several efficiencies desired from the rotor solution are found 





from; 
Т - 1 
μμ s Ás 227 
R Tee pn ee ( ) 
for rotor efficiency; апа: 
28 Δ 2 e Te, 
My аһ Te — This = 


for internal efficiency, The velocity coefficients are: 
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Computer Program 

The data reduction computer program, TTTRPGS, included in Appendix 
VI, proceeds through the various computations required in the solution 
of the foregoing equations by using subroutines described below, Inputs 
to the program are: the number of test data sets to be reduced; the date 
of test (month and day); ambient pressure and temperature; the number 
of data points for each data set; a value for the variable, TYPE, which 
indicates whether the data is for stator-only (TYPE = 0) or full turbine 
(TYPE Z 0) tests; and all experimental data, All experimental data inputs 
are the values as recorded during the test, except that temperatures are 
in °F rather than in millivolts, and rotor moments are in in-lb rather 
than in percent of maximum moment, 

Subroutines are called in the following order; 

1, FLORATE, Determines the flow rate measured by the flow nozzle 
(145) and the plenum labyrinth leak rate (170a) and (166). The turbine 
flow rate is then found as the former flow rate minus the later, 

2, 5ТАТОК, Solves the stator equations and determines the proper- 
ties at the stator exit. Values based on the solution of the momentum 
equation are computed and printed, but values used in subsequent calcu- 
lations are those obtained from iteration using continuity considerations, 
Also determined by this subroutine are the values for the area restriction 
factor, stator efficiency and stator velocity coefficient. Mean stream- 
line radius is taken as 4,26 in. 


3. ROTOR. This subroutine converts absolute flow conditions at 
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one to relative flow quantities, determines AH, and performs the iteration 
using continuity considerations to determine relative gas efflux angle 
and velocity, and temperature at the rotor exit. The other thermodynamic 
properties, flow angles, rotor efficiency and rotor velocity coefficient 
are computed also. A value of 4.45 in. is used for the mean streamline 
radius. 

4. REFERD. This subroutine converts data to their referred and/or 
non-dimensional forms, 

Outputs from TTTRPGS are the thermodynamic properties, velocities, 
flow angles for both the absolute and relative flows as appropriate, 


along with the non-dimensional values (k efficiencies, r*, etc, ), and 


ів? 
the actual and referred values of rotor moment, horsepower, RPM and flow 
rates, 
Results 

Three tests without the rotor installed were conducted on 27 December 
1965 and on 5 and 6 January 1966, Large fluctuations in axial force 
readings were noted during the first test and the data obtained were 
unusable, The turbine hood was removed to eliminate the backflow caused 
by impingement of discharge air on the hood, All subsequent testing, 
including tests with the rotor installed, was done with the hood removed. 
Data from the second stator-only test were not considered valid, since 
an open turbine plenum drain line was found at the conclusion of the test, 
Consequently, only the data from the third test form the basis for results 
presented for tests without the rotor installed. 

The turbine was rotated for the first time on 12 January 1966, and 
five additional tests were made between that date and 16 February. Rotor 
moments measured during the first two such test were obviously too high, 


since unbelievable turbine efficiencies resulted, Results of these two 
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tests were considered valid only for stator performance considerations. 
Turbine tests were discontinued during the fifth test run on 28 January, 
due to suspected bearing failure in the air dynamometer unit, Upon 
disassembly of that unit, the bearings were found to be satisfactory, 
However, the dynamometer rotor had been scraping against its housing. 
Moreover, the spline at the dynamometer end of the quill shaft connecting 
the rotor to the dynamometer unit was badly fretted, indicating possible 
misalignment. After reassembly and reinstallation of the dyanmometer, 
the turbine was realigned and testing was resumed. It was necessary to 
seat the dynamometer spline (male) as deeply as possible in the quill 
shaft to maintain adequate undamaged spline bearing and shear areas, 
Dynamometer bearing failure was again suspected on 18 February during 
preparations for the seventh test, At this time, the dynamometer was 
returned to the manufacturer for disassembly, inspection and necessary 
repair. Hence, testing was again interrupted pending return of the 
dynamometer and receipt of the new quill ‘shaft which had been ordered, 
The turbine became operable again on 15 April. However, it was too late 
for further turbine performance tests for inclusion in this study. 
Consequently, zero speed and zero torque tests were made as the final 
tests for consideration herein, 

It is pertinent to consideration of the test results that the total 
pressure probe pressure line parted inside the stator-plenum assembly at 
the conclusion of the third test, Values of Pt. used in data reduction 
for subsequent tests were determined by applying a correction, deduced 
from previous tests, to the static pressure readings, py, obtained at the 
probe location, 


Results of the stator-only and the turbine tests are shown in Figs. 
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4l through 49. For ease of reference, Table VII indicates the contents 
of these figures, Also shown in Figs. 41 through 49 are selected 
results of the performance analysis. 
Discussion and Conclusions 

As is evident from Fig. 41, stator-only tests yielded values for 
stator efficiency of 0.9, based on the continuity equation iteration, and 
0.93, based on the momentum equation. The difference between values of 
the gas efflux angle (OC) given by the two methods of solution (see 
Fig. 42) reflects the difference in efficiency values, since common 
values of tangential velocity were used in computations performed by the 
two methods, The values for eq used in the performance analysis are 
seen to definitely favor the continuity-based values, whereas the esti- 
mated stator efficiency used for the performance prediction was 0,9095, 
Any error in the continuity-based values can be due only to the stator 
exit area used in computations, This area was the stator exit annular 
area minus the product of blade heighth, number of blades, and the 
dimension for the projection of blade trailing edge thickness on the 
stator exit plane, Although accuracy of this exit plane area is dependent 
on manufacturing tolerances, the area must be considered quite accurate, 
The momentum-based values are very critical with respect to the pressure 
terms and also depend on the axial force. For a typical data point 
(рег /РІ ж 2,61), the fluctuation in the latter quantity over a period ої 
about ten seconds was $t 1,25 lbs from the average value of 24,65 lbs. 
With the flow rate of 3.44 lbm/sec, a variation of * 9.4 ft/sec results. 
Considering that the values of axial velocity obtained for this test 
point were 404 ft/sec from the momentum equation and 334 ft/sec from 
continuity considerations, the fluctuation in axial force measurements 


appears to be relatively unimportant. It must be considered then that 
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TABLE VII 


SUMMARY OF THE TURBINE TEST RESULTS 
PRESENTED IN FIGURES 41 THROUGH 49 


Figure | Shows | ὕετους | Remarks — 


curves for Pt,/P2= 
constant 


Referred Flow Rate, kis same as above 
/8/ 
Theoretical degree of Kis same as above 
reaction, r* 

















incidence | same as above 
angle, i 


curves for 
Pto/ P2 z constant 
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Stator Efficiency, Ns 
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either one or more of the pressures used in the momentum equation 
computations is in error or that one or more of the initial assumptions 
concerning the areas upon which the several pressures act is not valid, 
Subject to the limitation imposed by the assumption in the performance 
analysis of the loss coefficients being constant with radius, it can be 
considered that the pressure distribution from hub to tip given by the 
analysis would indicate that the average of hub and tip values of pressure 
should yield a representative value for the pressure іп the stator exit 
plane, pj. In addition, the pressure pattern given by the analysis--low 
at the hub and high at the tip--was observed during the tests, For 
instance, for a mean streamline stator exit Mach number of one, the 
performance analysis gave values (in psia) of 13.03 and 16.21 ас the hub 
and the tip, respectively, while matching test values аге 13,97 апа 15,84, 
Considering, however, that the performance analysis average pressure 
occurs at a smaller radius than that of the mean streamline, an average 
effective pressure obtained by integration of pressure times differential 
area would be higher than the average of the hub and tip pressure, It 
should be noted that such higher value would reduce the value of axial 
velocity obtained from the momentum equation, which is the effect required 
for stator-only tests, The entire stator exit plane annular area, 18,38 
in.?, is acted upon by the pressure pj. Fora change of 70 ft/sec in 
axial velocity, the difference between values given previously, a value 
of p] only about 0.42 psia greater than the average of hub and tip 
pressure would be required. Unfortunately, a pressure survey in the 
stator exit plane would be required to accurately determine the radial 
pressure gradient which would be needed for finding the force-average 


pressure, The performance analysis pressure distribution could be used 
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as an approximation, however, Further consideration of the pressure 
terms which occur in the momentum equation indicates that it would be 
desirable to know the pressure acting on the flow director бы! down- 
stream of the stator exit plane, on both the annular area at the sharp 
break at the stator exit plane and the annular area at the downstream 
end (see Fig. 39, p. 122). 

Any comparison of stator efficiency and gas efflux angle values 
given by the momentum equation and from continuity considerations for 
tests with the rotor installed would be subject to the considerations 
given above for the rotor. However, it is evident from Fig. 42 that for 
the rotor tests the values of axial velocity obtained from the momentum 
equation are too low (i.e., + too large). Of some interest is the 
fact that the radial gradient of pressure obtained with the rotor 
installed was somewhat more severe than that for stator-only tests. For 
the conditions considered in the preceding paragraph, the values of hub 
and tip pressures measured during tests with the rotor were 13.5 and 
16.41 psia, respectively. These considerations serve to emphasize 
previously-made comments that the pressure acting on the cover plate 
must be known, 

The parabolic nature of the equation for the radial distribution of 
pressure acting on the cover plate indicates that pressure taps should 
be installed at the center of the cover plate and at a radius of 0.707 
times the cover plate outer radius. The redesigned thrust and torque 
flexure will, of course, yield the force acting on the cover plate, as 
previously discussed. Even though this flexure and pressure taps in the 
cover plate effectively give the same information, installation of both 
is considered desirable. 


It should be noted that the spread of stator efficiency values for 
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tests with the rotor installed, Fig. 4l, is not reflected in values of 
C, Fig. 42. This situation is considered to be caused by the difficulty 
noted previously of obtaining all readings rapidly enough to preclude 
effects of variation of test conditions. Since provisions for automatic 
data logging are nearing completion at the time of this writing, this 
problem should be alleviated in the near future, 

Of considerable interest is the increased stator efficiency with 
the rotor installed, which can be estimated as 0.93 to 0.94,over that 
without the rotor, This result belies the assumption that conditions 
downstream do not affect upstream conditions. Since it seems unlikely 
that stator blade profile losses are reduced by the presence of the rotor, 
it is conjectured that the rotor acts to reduce the effects of secondary 
flow losses, Considering the sharp break in the gas path in the stator 
exit plane, it seems probable that the presence of the rotor may provide 
sufficient flow guidance in that area to reduce losses due to vortex 
formation as the flow spills over into this relatively large region, 

The increased severity of pressure gradient in the stator exit plane for 
tests with the rotor installed over that observed for stator-only tests 
indicates that the flow pattern in the area is altered by the presence 
of the rotor. Installation of pressure taps in the flow director ring 
at this break could, perhaps, give some insight into the nature of the 
flow alteration, as well as being desirable for application of the 
momentum equation, 

The curves of Figs. 43-46 illustrate the performance of the turbine 
as measured in tests three through six with the rotor installed and with 
an axial clearance of 0.243 in, Some evidence of the effects of non- 
simultaneous readings exists in the scatter in values of internal effi- 


ciency, Fig. 43, whereas the moment curves of Fig. 46, which depend only 
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on flow rate and measured torque, do not exhibit undue scatter. Constancy 
of referred flow rate is indicated by Fig. 44, with wpef 5 1.375 lbm/sec, 
The variation of theoretical degree of reaction indicated im Rig. 45 is 
considered quite reasonable for an impulse-type turbine. Figs. 48 and 

49 are included here for illustrative purposes and discussion in the 

next section, with respect to comparison of performance prediction and 
test results, 

Both internal efficiency, Fig. 43, and rotor efficiency, Fig. 47, 
vary with incidence, angle in a manner characteristic of an impulse-type 
turbine, i.e., with a narrow range of incidence angles for high effi- 
ciencies, Except for the differences between maximum rotor efficiency 
values for the different pressure ratios, the curves of Fig. 47 could 
be replaced by a single curve, 

The values of rotor gas efflux angle are not depicted graphically, 
since this angle varied over a small ranges -67,1 to -67.8 , with the 
higher value applicable to higher pressure ratios, Pc / P2» In addition, 
the stator area restriction factor, Е » was essentially constant, vary- 
ing from 0.944.to 0,952, and with the variations indicative of no dependence 
on stator pressure ratio. 


7, Discussion and Recommendations. 





Comparison of experimental and performance prediction results 
indicates that the most serious deficiency in the performance analysis 
is in the variation of the area restriction factor, Fig. 20, while 
experiment yielded an essentially constant value, The effect is indicated 
in the reduced referred mass flow rates after the flow is choked given 
by the performance analysis, Fig. 19. The effect is evident also in the 
agreement between the predicted curve of referred axial velocity and 


experimental values at lower pressure ratios, Ρερ/Ρι» with subsequent 
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departure at higher pressure ratios, Fig. 48, Apparently Ё should be 

taken as constant after the flow is choked, so that W ef Will be constant, 
Quite obviously, any performance prediction is no more accurate 

than the values used for stator and rotor efficiencies, The differences 

in magnitude between predicted and observed values of internal efficiency, 

rotor efficiency, and stator efficiency with the rotor installed serve to 

emphasize this fact. The high deflection through the rotor, 129 , makes 

estimation of Te quite arbitrary, since very little is known about 

turbine losses for rotors with deflections of greater than about 110- 

115°, General practice in gas turbine design limits the turning angle 

to 105°, Predicting the positive stalling incidence angle is subject 

to the same difficulties attending the choice of rotor efficiency, From 

Fig. 47, the stalling incidence angle for the transonic turbine can be 

estimated as about five or six degrees, whereas the predicted value 

used was 13.4°. The effect of this difference is evident іп Figs. 43, 

47 and 45. The much greater low-loss ranges for the prediction curves 

of Figs. 43 and 47 are clear indications, However, the manner in which 

the theoretical degree of reaction is influenced (Fig. 45) requires 

further explanation. One must first imagine that the predicted curve 

of Fig. 45 is moved upward so that its position corresponds to the higher 

efficiency (over predicted) obtained from the tests. Then, assuming 

that the loss coefficient remains constant, whem positive stalling 

incidence is reached, the smaller the value of this incidence angle, 

the higher will be the RPM at which the loss coefficient becomes constant, 

This consideration means that the negative reaction portions of the 

prediction curves of Figs, 29 and 47 would flatten out (towards zero 

reaction), 


The gas efflux angles used in the performance analysis proved to 
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agree very well with those determined. experimentally, Fig. 42 indicates 
almost exact agreement for the stator, while the predicted values for 

the rotor angle of 67.76° to 69 agree favorably with the values of 67.1 
to 67.8? observed experimentally. Further, the exit Mach Number for 
which the value of 67.8 was obtained was about 0.6, and the corresponding 
value о/з used in the performance prediction was 67.85^, 

The discrepancy between predicted and measured values of referred 
rotor exit axial velocity, Fig. 49, can be attributed largely to the 
differences noted previously between predicted and experimentally 
determined values of the stator area restriction factor, 

It is considered that the addition of instrumentation recommended 
earlier in this paper is necessary prior, to further testing. In 
addition, it is recommended that attempts to correlate the performance 
prediction program as given herein (or a similar program) be continued, 
by assuming different loss coefficients and different radial variation 
of such loss coefficients, It is pertinent in this regard that a thesis 
reporting on cascade tests of the converging guide vanes and the double 


circular arc rotor blades will be published shortly. 
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APPENDIX I 


FORMULA DEVELOPMENT FOR TURBINE PERFORMANCE PREDICTION 


А» Equation of Motion for Relative Flows. 
| The equation of motion for relative flow in vector form, (2), 


| р. 39, 15: 
VHa- Wx(VxW * zw) * TVs (2) 


Expanding the several terms of (2) gives: 





























roo ' *oz ar (231) 
W -- [ lp | _ | 0M. DW, 
ve) les dE. то Зою я? 
9 r 9. | асо " 95η T 
БТ: Oz or y \or 96 
М ni А κ. 
ет Wa Nr гл д2 
(232) 
Wx(VxW)— le ГА tr 
Wu Wa Wr 
(3 - 9) (awa ам) 12717127 
д ғ ar r\ or 90 r oo FE: | 
T i[owa  o(rW O(rWu) | ow, 
ове свечу) 
дум o Wa н дм δίτ Wa) 
ЕСЕН ο w + 9 - см) + 
| | a (r WA) oW oW Ə Wa. 
r | wee p 20 | - Μα E ст ] (233) 
Wx2ü z(igN,* 3M, +i, W.) x (ia 207) = le (2054) — Le(2 ax W,) (234) 
T Vs = f (ier $$ 131 95 +i OS ) (235) 
i fr 968 т ог f or 
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Combining (231) through (235) and writing the equations for each unit 


vector? 


ы ы» - йа), СШ, м, "НТ. (амь 


Wu E ~ 2 Wal 4 т 23 (236b) 


(с -ӘШЕР жи. [SCRI i x, [9^ -. 3948]. zor W. +T 35 (236c) 








For assumed axisymmetric flow, all derivatives with respect to Oare 


zero. Then: 


l: O=— We Ə(r W.) _ Ww. olr WwW) 
































ұу Т ar 2 a WW (237a) 
|. дНа _ Wr 2 W, шей Wo) E 
i: S = We ole) ος ὅδε 2 uw Т 25 (237c) 
Eq. (237a) can be rewritten: 
Ə(rW.) _ _ wç Ə(r W.) _ Wr 
2% wor τω (238) 
Substituting this expression in (237b): 
2 Ha _ We κ нь - We δ (τ ΜΙ Wa Wr 95 
ὃν = Мда Νε Man ЭР - Wa me 
ог: | 
ӘН» __ _ Wa W.W. an ) o Wr JW: " Wa We 05 
да — rW, or σαν ror MEE Wa. pu: БЕ мэн 


Multiplying (237c) by W, and (239) by Wa: 











H W. r (r W.) κ OW, oW, 
W 9 Ha C Wu We e OLEWA) wy er — WW, a Wa τὸ 
^92 | Fr or 4 дг 93 
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Adding these two equations: 


м, ЭНЕ +W. eB = T (w, 98.4 4,95) . (240) 


For adiabatic flow, total relative enthalpy is constant along any given 


streamline. Hence, the stream derivative is zero. With Ə/99 = 0: 





— 3 Sc H H 

МАЕ УН, = О E W, 2-5 + "E | (241) 
ог 

ota — - - (2414) 


252 М205 (242) 


Substituting (241a) and (242) in (239): 


ме ӘН мМ (ИМ) 9We. w 9Wa ; ΜΜ τὸς 
Wer rw, or Wos “esr ^9 we Ww. ar 


Then, multiplying this expression by (-W,/W,) gives: 











дНа _ Wu (rf Wu) дү 2 W 2s 
warar М iW QS 5 (243) 


Comparison of (243) and (237c) shows the expressions to be identical, 


Hence, this expression is the equation which must be solved, Rewriting 


fat 
(243) and with (м, 9% 120%) there is? 


м u ы 
8042) 2 зүү dW 4 2Wa MWe). 4 oW 2 Нв anaE o (244) 


Eq. (4), p. 40, can be written: 
2 


H. = He — = (4) 
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Also, (5),. p. 40, can be written: 


2 


Ё 


mE. m 
D 


(9) 


һә 


Ср 
Eqs. (4) and (5 ) are now introduced into (244). The nature of the 
equations relating total relative and equivalent enthalpies dictates 


that (244) be applied to the rotor discharge with this substitution. 


Then: 
οίννα͵) _ 3 W- м, ӘМ.) а. οι у? 
ασ δω. 235 (Pe - 2] 
H, _ ΜῚ 95 
+ 2/ esee (245) 
The multiplier of the final term can be expressed: 
| 2 | 2 1 We, 
2 1) | " б А 
Cp 2 e (Wa, + Wu; + Wr, = Ср |24. Wa, Wa, ( T Me) (246) 


Eq. (246) becomes, with (tan*A = we /w2) and (l+tan* A =1/cos*,_ ) 


Cp соѕ?А, 


рн. == „гк W, ж... (247) 


Ср 2 


Eq. (245), combined with (247), and with some rearrangement, gives: 


2) 2 
oM Ул 7 үү ӘМ, E Μία; 95, |. 2 Wa; m 44 «Ма, 





4) Ὁ 2 Cp C0S*A, OF r; 
aH. , a(Ui), 1 ЕР, M 


There is also: 


2 Er: 
ου). доме) M 2 олт, 90. ERON (249) 
2 2 





Еа. (248) and (249) combine to give: 
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w r 2 (3427 
ХӨ up μμ й 


a JZ ερεο52λ;, OG νου 


-2ата2|2ң,-М93-0 (250) 


To make this equation non-dimensional, it is multiplied by (ενα), 


where the second subscript, m, denotes the mean streamline value, Then: 


2 W. 2 
5 9 (м) = 2 Wa, Waz A e (Ws) edi Way Cam 95, 159 WaW, m ЭС. заг ον) 
Warm dh М Wa, *m дж М Cp C0S7A, OV, Waren Λα οί см) 





Oo Wa Wa = ба ӘН, 2 от 11 Їл [zm fam 2 He = μιὰ Wa, OS: 251 
"a Wa, . Wa, ^ WA Or б LA M Ср А И Was 9f, á | 


Here, it will be observed, the subscripts, 72!, can not be neglected un- 


less the equation is to be applied to that station (rotor exit). Since 


the station of interest is station two, the subscript will be dropped. 





Let: 
- _We 
Y WIL (252) 
Х- = (253) 
and: 
6#— -ᾱ- (254) 
(p 


Eq. (251) can now be written in terms of these non-dimensional quan- 


tities: 


EN ; 


ox ¢ 





y 9M γ΄ Os an oY хапа ΟΙΧῪ Μα) 


2 соз?х ӘХ ОХ 


El 
QJ 


Α Um Y tan/ó _ 2 ah ELA Um Uz [ο -Y'tajg t - O 


Wam Wy, OX Was, — |Wa, (255) 


But; 
tans t Y X 
2Y P T zy tana xy den хавгүй +Ytana 2X 
2 2 
= “tana Stang + 7 Ytanzg SY + 27 Сап,” 


Also, there is: 





Atang _ _! af 
aX C0824 ӘХ 


апа: 





Еа. (255) then is: 


Y W з 95" 2.4 д EE 
3 20 + ап) - -2у r 1-9 +2ү 2412.96 + i 





ин Ytan® . E ОН, аа у" Кап (95 (256) 


ат 


Multiplying by (cos ^9 /1*), and, since (1+ сап 7 = 1/соз%4 ): 





Wa OF -imat 2tany A +$ sing + 


E μμ ο cos _ 2016 ӘБ: + EHE CO sinala =0 * 49257) 
EY Wam Y ðX Wam Y p^ 


The second term of this equation accounts for the effects of streamline 
2 
curvature, as represented by the cosZA and (- = 2 W ). Tig. 15, 
a o2 
p. 42, shows the sign convention for streamline curvature (i.e., positive 
for the radius of curvature in the positive coordinate direction). Fig. 
15 defines also the streamline displacement, r, and the length over 


which this displacement occurs, L. Streamline curvature can be approxi- 
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mated; 7 


б В янь x SE ; K&4 tob | (258) 
2 m Р 


It must be noted that (258) and (259) are defined here for curvature 
across the rotor only, because of the small blade heighth of the stator 
and because of the large flare area available across the rotor, In 
(258), k is positive for positive curvature. Putting (258) and (259) in 
(257): 


diny’ * d- (+ 24 rp &p) (Es) de _ 2tang JÉ -Z sing — 


ο ο eo 


It is evident that (260) is non-dimensional, except for the terms contain- 
ing equivalent enthalpy in BTU/lbm; let: 
СЕ 265 Т ( ££?-1bm/BTU-sec?) (261) 
Then, (260) becomes completely non-dimensional as: 
а — 2 Ч 2501 - (£2 ds*| _ ід 2 сіна 
jx = — С052)|-| ХК ү Ъ |9 Ztane É X 5112 — 
40.5112 со59 _ 2010 со522 PEDES fate coge Ci He cos’? _ sin vel ds" (262) 


Wam Y м ү? Wan ¥* а м2 ү? 
Κα. (262). the final result, is included in the basic paper as (26). 


17 үдуга, M. H. Aero Thermodynamics and Flow in Turbomachines 
(New York, London: John Wiley and Sons, Inc., 1960) p. 453. 
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18 
В, Theoretical Basis for the Area Restriction Factor, 6 : 


Fig. 50 shows typical conditions 
at the exit of a blade row, with boun- 
dary layer on both sides of the flow 
channel, Outside the boundary layer 
the velocity is very nearly the the- 
oretical velocity. Within the boun- 


dary layer, turbulent flow can be as- 





sumed, so that the velocity ratio can 

Illustration of Boundary 

be taken as: Layer Effects at the Exit 
of a Blade Row 


ως 2 26 F 
v= (4) р та (263) IGURE 50 


The mass flow rate leaving the blade row can be written: 


5 
M = Aru Vru 205 Ха [s- Exxy — PT uye tiy (264) 


Here, fru and Vry are the values for isentropic expansion from inlet to- 
tal pressure, Pto» to the pressure at the discharge, pj, which is assumed 
constant across the spacing, s. 
Fig. 5l depicts the thermodynamic 
process for the flow outside the 


boundary layer, and for a point in- 





side the boundary layer. It is evi- 
Thermodynamic Process 
dent that the first term of (264) through a Blade Row 
gives the flow rate outside the boun- FIGURE 51 
dary layer, while the sum of integrals is the flow rate within the boun- 
18yavra, M. H. Problems of Fluid Mechanics in Radial Turbomachines, 


Parts III and IV (Rhode-Saint-Genese, Belgium: Von Karman Institute for 
Fluid Dynamics, 1965) VKI Course Note 55b, pp. 646-50. 
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dary layers. The discharge angle is given by: 





Е а 
ба COS з тент) (265) 


With (265), (264) becomes; 





x 

T acosa] [scosas -t 5 созой =) Lud 

m = fra Ven coeds als Cose - #5@ = V ZA Wa 4 ra d 54 ун Мен 4 
Cancelling terms, there is: 


rn = Aw a| i= $ À (- [Z щі) (266) 


For the constant pressure assumed at the discharge: 


ο ων τι μα. _ _1\—Хе (267) 
T Т. (п — Ten) (LY | — Xe(& Y 
УЗ to (т, ) (1) 4ГА 
where: 
Y-1 


es (317 (268) 


Eq. (266) can then be written, with (263), as: 


үл = д, ma f1- Е (1 -X0 [2 sen] (269) 


Defining the so-called displacement thickness: 
+ | n 
"= 5|1- (Ёл Ke) | 25у dn η (270) 
Eq. (269) becomes: 
4 
M= ANu A [ -- zs (271) 


For Уа, the actual average discharge velocity, the loss coefficient 


of the blade row is: 


m (ay 
= σα (272) 
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This loss coefficient can also be expressed in terms of the average loss 


of kinetic energy. Hence: 


DE. a s "ЗЕ 273 
Чийгээ чийг Ч -- 


where E is the actual kinetic energy of the flow: 


2 $ i 
E=2,V,,(a-< 5) + S [ оо e dy (274) 


Again, the first term represents the value for the extra-boundary layer 
flow and the second term the value for flow in the boundary layers, 


Eq, (274), transformed in the manner used to obtain (269), is: 
un 5 he 
E= Ay а-о “уяач (275) 
Then, the so-called energy thickness is: 
Хы 
д h- а-ә 6 Ж ағ dn | (276) 


The loss coefficient can now be expressed, using (271), (275), and (276), 


as; 






PUT ll s i d 


Arn Sea I-32 





ог: 


5 
ееш шинийн (277) 
| -2 а 
The area restriction factor, Е » represents the fraction of the flow 


area which would be usable for uniform theoretical velocity. The defi- 


nition of € is then: 


É = | — 55 (278) 
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Also, the energy parameter is defined as: 


н“ BL | (219) 





Using (277) and (279): 


444 
554 


E H 


Ie o e s s. 280 
RN — ee “2 


Неге, f is the profile loss coefficient. 


| 


C, Development of the Energy Parameter. 
The integrands of (270) and (276) can be evaluated by expanding 


the denominator by means of the binomial theorem: 
=I T m 3 
(1 Xe") що exer Xen + Xe” + XEN... 


The integral of (276) is then: 


СЭ ТЭНХЭЭ +. dn 








Integrating: 
і MS | Χο Xe Ke Die 
| l- Xen?” шэн jm*l Ш бүү? | Ч m "inl" eo 


Еа. (276) can now be expanded, giving: 











чаї 4 2 3 4 5 
Or = |-(-.- + qh + i, Hy хао оха ae, 
$ m Эт! біті Imel fei Пи + і Зт+ | бт + Ї Төзі Ямы © tims 
ог; 
ке 2 3 4 
б = ГАК) ЫБ НЕБЕ κος Ae) 
5 Зт +! Smt! 7т +! Өт + |! ll mti 
Then: 
РҮҮ, 2 3 4 
6 _ | | Хе Хе Хе Хе ) 281 
5 2 (Хе-1) Υρ--| T 3m+l x 5m +l t Tm + x От + | i Ше +1 ( ) 
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* 
By similar expansion and integration, there is for 52 





— + 


- 2 3 4 
| | 1 
— mei Xe h- -0 в (282) 


OR a 2 
о |) сто о, т. е 


Then, from (281) and (282), the energy parameter is: 








Xe + Xe + Ж E h 

















| | 
ае 1 с> ЙА... 
44% -| 4 
y- Xe 3m*l 5m +! 7 т+! 9 т! ит +1 (283) 
| | * 2 A о 9 , 28 
Хь – | m t I 3m | 5m *l 7m | 9mti 


Fig. 52 shows H*** as a function of pressure ratio;for Y - 1.4 and m - 0.15. 
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APPENDIX II 
COMPUTATION OF FLOW ANGLES AND LOSS COEFFICIENTS 
A. General, 

1. Тһе basic method used to estimate gas outlet angles and loss 
coefficients is based on formulas and graphs developed by Ainley and 
Mathieson. [4,5] Where computations are based on information from other 
sources, footnotes will so indicate, Fig. 53 shows definitive blade 
geometry pertinent to the calculations. 
B. Formulas Used. 

The General formula for gas out- 
let angle prediction for Mach Numbers 


of 0,5 or less is: 





o 3 tan'[(-4 tan oq + Атап 21 (284) 
Where: 
BLADE GEOMETRY USED IN 
GAS ANGLE PREDICTIONS 
A=F A cost. сөзді. (285) 


FIGURE 53 


Here: F is a factor which accounts for the type clearance (1.35 for the 
unshrouded blades of the TTTR rotor); h is the blade heighth (in.); Οἵ. 
is the inlet blade angle; and X is the gas outlet angle, uncorrected 


for tip clearance, given by: 
* 5 
Οία = X, t 4$ (286) 


In (286), s is the blade spacing and e is the curvature of the back of 


the blade; These two quantities are given by (see Fig. 53): 


DC cos" | E (287) 
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and: 


ie | 
е = b (288) 





When a blade row has zero tip clearance the gas outlet angle is given 
by (286). 
For a Mach Number of one, the mean streamline gas outlet angle is 


determined from: 





ο), = cos M (289) 


where Ав, is the throat (exit opening) area and: 


nn) SS s Up -D*2K (290) 


It is evident that (289) is basically (287). Flow angles for other Mach 
Numbers (between 0.5 and 1.0) are estimated by drawing a smooth curve 
between the values for 0.5 and 1.0, but with an inflection point at 
M = 0.75. 

The basic equation for the profile pressure loss coefficient for 


zero incidence, (Y5)i:.o, is: 


- 


27 
= 4 ὁ (x) [ а и а]! a (291) 


Here; t is the maximum blade thickness, and the ratio, t/c, is restricted 
to a maximum of 0,25; and the equality subscripts indicate which data 
curves are to be used. The values of profile pressure loss coefficient 
for other then zero incidence angle are obtained from зару wiso 
versus i/i,, where i/i, is the ratio of incidence angle to positive 
stalling incidence angle. This latter value is defined as the positive 
incidence angle at which the profile loss is two times its minimum value, 


For the calculations, i, is determined from curves based on a solidity of 


S 
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0.75, with corrections applied for the actual solidity. 


Secondary and tip clearance losses are determined using: 


\„+\,=[х Ч | хози. 





(292) 


where: CX , the mean of the gas inlet and outlet angles, is a function of 


incidence angle; B is а type-clearance correction factor (0.5 for radial 
clearance); and: 


T Е (tan αι * tan e) cos o, 7 f(i) 





(293) 
X z £[(A/A) | (109 1070.07] (294) 
Here; 
АЕ ра (295) 


The terms, I.D. and 0.D., are the annulus inner and outer diameters in 
the exit plane, respectively. Values of ^ are obtained from a curve 
based on the bracket quantity of (294). 


Leaving losses are found from: 19 


2 
А сағ (8/5) 


(1 878) Жэ 
where, see Fig. 50, p. 156: 
9 te 
з — — 
ee (зт) 
е 59 
апі: 
sS “a мі (298) 
a Н” э | 3- ($,* $, « 5.) 


l9?yavra, M. H, Aero-Thermodynamics and Flow in Turbomachines 
(New York, London: John Wiley and Sons, Inc., 1960) p. 89. 
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The energy parameter, НЭЭЖ, is obtained from Fig. 52, p. 161. 
The pressure loss coefficients given by (291), (292) and (296) are 
defined as: 


„ Loss in Total Pressure 
Outlet Total Pressure - Outlet Static Pressure 


For use in the performance prediction program, values obtained must be 
converted to the loss coefficient, Р . For incompressible flow, assumed 


for these calculations of loss coefficients, the required conversion is: 





ү (299) 


This conversion should not be made, however, until the total pressure 


loss coefficient has been determined, since: 


xt; n Ἔ Tov (300) 


Nevertheless, in these calculations, it was necessary to convert to £ for 
the sum of profile and secondary and tip clearance values of Y in order 
to obtain the leaving losses. No appreciable error results due to the 
relatively small magnitude of such leaving losses. 

С. Calculations. 

Calculations for determining stator gas outlet angles are given in 
Table VIII, while stator loss coefficients are computed below, Table IX 
contains the computations for the rotor gas outlet angles and those for 
the rotor loss coefficients are tabulated in Table X. 

Figs. 54 and 55 show the variations of predicted stator gas efflux 
angles with Mach Number and Radius, respectively; Figs. 56 and 57 show 
similar information for predicted rotor gas outlet angles. Fig. 58 shows 
rotor total loss and velocity coefficients versus incidence angle ratio, 


і/і се 
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Estimation of Stator Loss Coefficients 
Losses will be taken as constant with radius and equal to those for 
the mean radius. 
ЗУ ой, ой): 
6/с-0,603; О(-0; OX, =76.26° 
Ы ий. (from curves) 
Y. 7 A(C, / (s/c))" (cos^o« /cos? ox. ) 
και etan 1/2) =63.95° 
С, / (s/c)=2tanc&, со5 =2(4.09) (0.439) =3. 59 
(с./ (в/с))4512.9 
=# (А, /A,)*/ (141.5./0. D. )) 
A, ,7aZh=(0.197) (31) (0.69) =4.216in.” 


2.7, 79° ae. Jen. - 


A, 7(TI/4) (9.17 
(A, /,) ^ / (1H. D, /0, D, ) (4. 216/18.38)? (147. 79/9, 17) -0.02855 
-0.0054; cos0(-0.0557; сов30<, -0.0848 
Y (0.0054) (12.9) (0.0567/0.0848) -0.0467 
Кы 
$ +$ =(0.0967/1.0967)=0.0883 


À =cos“o(, ( Š iym. as 


mix 
н “5.2.28 (р/р, =0.843); t =0.024 
6 a C5, $ rai ae $,* 5 )=(0.0883)/ (1.2840, 0883) =0. 065 
6/5=( 5 /a*t /a)/ (1 + t /2)-(0.06540.1218)/ (1.1218) 0. 1663 
A 1,7 (0.0563) (0.1663)7/ (0.8337)? - 0.00225 
f = 0.002245 


Б окат ^ 0.09054 
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TABLE VIII 


CALOULATION OF PREDICTED STATOR GAS OUTLET ANGLES 





radius (r) - ín. 

radius ratio, X 2 г/г, 
chord (c) = in, 

spacing (s) - іп, 

exit opening (a) = in. 
blockage (5) - in, 

ez 32 /8b; bz0.16; j«0.76 
8-5, 


сов”? uou) eades,. 


зо 
«Үү = ?(сов”1 23) 


“ti, 
4(s/e) 
“(үк «ү + (а/е) 


нв. MEAN ES 
М < 0.5 

5.895 4.24 4.585 

0.919 1.0 1.08 


uu oer rs THIS ша. 


0.7895 0.8594 0.9293 
0.177 0.197 0.223 
ын жез pop ea a 
"------.-.... 4,52 ----------- -- 
0.7185 0,7684 0,8585 
75.75 75.50 74.95 
75.73 75.50 74.93 
0.698 0.76 0,822 
N= 1.0 


Ath = (31)(0.197) (0.69) = 4.216 in? 
À, =(m/h)(9.172 ~ 7.792) ~ (0.071)(51)(0.69) = 16.86 %а2 


— 51 


Cu s 1 


= соз"! (4,,/4) - сов” (3.216/16,86) 


15.5 


AX: (Οἵ; )м E] = (с, )κ«ο.6 - -0.76 


Οἵ μα 1 
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Stator Gas Outlet Angle -deqg. 


Stator Gas Outlet Angle -deq. 
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VARIATION OF PREDICTED STATOR 
GAS OUTLET ANGLE WITH MACH NUMBER 
FIGURE 54 


4.0 4.2 4.4 
Radius-r,-in. 


VARIATION OF PREDICTED STATOR 
GAS OUTLET ANGLE WITH RADIUS 


FIGURE 55 
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TABLE IX 


CALCULATION OF PREDICTED ROTOR GAS OUTLET ANGLES 


ITEM 
radius (r) - in. 
radius ratio, X = r/r, 
chord (c) 245. 
spacing (s) - in. 

exit opening (8) - in. 
blockage {t+,) - іп. 
(s-t, ) 

сов"! (шаа... 


) - дек. 

Be z ο (----2---)) «Өр: 
Я Sete, 

4ο = 45 (since 4s/e = 0) 

tip olearance (k) - in. 

blade heighth (h) - in. 

k/h 


type-clearance factor (F) 
cos ó; (£ = 62°) 

сов, 

tan 22 

À = P(k/h)(cos #)/cos 5 ϱ) 
(1-4) | 
(1-A)tan Z; 

A(tandi) 

O = ((1-A)ten Z 2 * А(їап,®,)) 
6, = tan”? (0) 
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HUB MEAN ТІР 
5,825 1,55 4.875 
0.88 1.0 1.12 
-—— τ е=— 
0.4006 0.4555 0.5224 
0.117 0.157 0.166 
ιο о ка 
0.5409 0.5958 0.4627 

. «69.95 -69.7 “69.0 

-69.45 -69.2 -08.5 
-69.45 -69,2 -68.5 
---------- ---0.025-------------- 
-------------1.05 -------------- 
-------------0.0258---------.... 
-------------1.55 -------------- 
-тесес------ -0.469 ------------- 
0.551 0.255 0.266 
-2.67 -2.61 -2.54 
0.045 0.0425 0.04125 
0.957 0.9575 0.95875 
-2.555 -2,525 -2.455 
0.0809 0.0799 0.0775 
-2.474 2 -2.445 -2.358 
-68.0 «617.16. 291.0. 


- шел ы "—— “шин 


TABLE IX (continued) 


ITEM | HUB MEAN 


The foregoing calculations are for М < 0.5. 
| zl, Man Radius 
A, = A (l-k/h) + Dk 
aZh = (0.137)(60)(1.05) = 9.45 
(9045) (0.9762) + (9.75)(0.025) = 8.43+0.0766 
Aq 5 9.196 in* 


> 
u 


М 


Аъ, 201/8) (О, + 2к)?-0б) - ++„2һ 

=( 1/4) ((9.8)? = (7.65)*) = (0.0579) (60) (1.05) 
Ano = 29.2 = 5.645 = 25.655 іш? 
So = сов”1 287729, = 00871 (9.196/25.655) 


/З 2 - deg. -69.0 


4,4 2 = AP за AP 505° -1.24 deg. 


ΑΦ - deg. -09,24 -69.0 
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~58, 24 


іт 2 -(deg. ) 


Rotor Gas Outlet A 





Mach Number 


VARIATION OF PREDICTED ROTOR 
GAS OUTLET ANGLE WITH MACH NUMBER 


FIGURE 56 





Rotor Gas Outlet Angle, 4$? -(deg.) 


Radius, rj.-(in.) 


VARIATION OF PREDICTED ROTOR 
GAS OUTLET ANGLE WITH RADIUS 


FIGURE 57 
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TABLE X 
CALCULATION OF PREDICTED ROTOR LOSS COEFFICIENTS 


Profile Loss Coefficient, Zero Incidence. i A, ^ 
| ж 
а ' 2 ; t/o 
(а= oy [opa 0) * 2) M =) ` "виз, - NIS 
{/ο 5 0.296/0.751 = 0.394 use 0.25 
/ а 629 
/2 2 -69.2 
(21/62) = 0.803 
з/с = 0.4555/0.751 = 0.607 
Ү,, үз 0.14 
473) 0.897 


1»(1- 0) z (0.0584 0.805(0.14 - 0.038)) (0.25/0.2) 


b. Stalling Incidence (i,) 
owe - z 1.05 
/ Se (s/o « 0.75) 
o Sa = -—57.76° 


Alajo.) ^ "9*5 


415 = 1. - (в/с = 0.75) 5 4. 
/ VI (s/c » 0.75) * 70-982 


4 А. z 9.09 (value is conjectural due to extra- 
8 (s/c —0.75) polation) 
i am 20247 
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TABLE X (continued) 
ὃν Secondary and Tip Clearance Losses. ! 
TEX =[a +B(k/n)] [ 0,/(8/c У12 (сов2,6,/сов24,) 
В 20.5 
k/h Z 0.0258 
cos d, z 0,1435 
АЕ (а/м Уа T I.D./0.D.)) 
A Z Zhao = (60)(1.05) (0.137) * 8.65 in^ 
А12 Zhya, = (60)(1.05)(0.194) = 8.85 in? 
(1 1.D./0.D.) = 1 +7.65/9.80 = 1.165 
(а/м ^ /(1--1.D./0.D.)) z (8.62/8.85)- (1/1.765)2 0.538 
A= 0.03 (extrapolated) 
Bk/h = 0.0119 | 
A+Bk/h 2 0.0419 
(A+ Bk/h)cos?4 5 x 0.00601 
Y,- Y, s 0.00601 (01/(а/с))“(1/сог22 |) 


1 -40,2 -20.1 -10.12 0 6.1 15.4 
tan £ 0.4 0.856 1.247 1.88 2.565 2.84 
(tend, * ten 45)/2 -1.022 -0.7875 -0.599 -0.2825 0.6 0.1975 
din δ «διό  -38.5  -30.95 -15.8 5.45 11.2 
cos Z 0.7 0.785 0.857 0.962 0.997 0,98 
tan J, = tan ó; 2.845 5.511 3.692 4,525 5.010 5.285 
о1/(8/о) 2.98 5.20 6.52 8.55 10.0 10.54 
8,/(8/с))2 15.9 270 4.0 69.5 100. 107.5 
(1/вов2/4,)) 2,92 2.07 1.588 1.122 1.008 1.065 
Y + Y, 0.279 0.338 0.382 0.468 0.606 0.689 
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i/ia 





T Yyy 


TABLE X (continued) 


à. Total Loss Coefficient. 


ἵρ/ (i =0) 


Y 
P 


Y TY ,+ I 


bot bat Sy 
-ἲ (50.2) 


вам 


Β 


жи 


5"/а 


t/a 


6 /a + t/a 


S/s 


(1- 5/9) 


(5/4) 
(1- $/s 


( S/s)*/(1- 5/8)? 


A mix 
b nix 


$ total 


n=l- 052 


Y=V4 


Sp 


j^ 


0.279 0.338 0.382 0.468 0.606 0.689 
4.28 2.1 142 1.0 1.22 2.08 
0.626 0.3075 0.208 0.1464 0.1787 0.3045 
0.905 0.6455 0.59 0.6144 0.7847 0.9935 
0.475 0.3925 0.371 0.38 0.439 0.498 
1.28 — 1188 1.26 1.8 1.4 1,28 
1,759" 120729 1.651 125 1.71% 1.778 
0.271 0.235 0.225 0.229 0.255 0.280 
0.1824 -----ваше------------------------.--... 
0.4534 0.4174 0.4074 0.4114 0.4374 0.4624 
0.5855 0.353 0.345 0.348 0.370 0.3915 
0.6165 0.647 0.655 0.652 0.630 0.6085 
0.147 0.1245 0.119 0.121 0.137 0.153 
0.58 0.419 0.429 0.425 0.396 0.57 
0.587 0.2975 0.278 0.285 0.346 0.414 
0.0555 0.0426 0.0399 0.0409 0.0496 0.0594 
0.0526 0.0408 0.0384 0.0393 0.0473 0.056 
0.5276 0.5555 0.4094 0.4193 0.4865 0.554 
0.4724 0.5667 0.5906 0.5807 0.5137 0.446 
0.688 0.752 0.77 09046: 0.717 0669 
Take the profile loss coefficient as: 

ὃς (αρ ΕΣ ту) Protas 
0.5295 0.187 0.1509 0.0906 0.10 0.1529 
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Loss/Velocity Coefficient 





Incidence Angle Ratio - i/i, 


VARIATION ОҒ ROTOR LOSS AND VELOCITY COEFFICIENTS 
WITH THE RATIO OF 
INCIDENCE ANGLE TO POSITIVE STALLING INCIDENCE ANGLE 


FIGURE 58 
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APPENDIX III 
COMPUTER PROGRAM FOR TURBINE PERFORMANCE PREDICTION 

A, Inputs. 

Listed below are the inputs required for Program PERFORM in the order 
of the READ statements. 

1, Stator mean streamline gas outlet angles (ALFAM--radians) for 
stator mean streamline exit Mach Numbers of; 0,5, 0,7, 0,75, 0.8, апа 1, 

2. Same as (1.), but for the rotor (ВЕТАМ), 

3, Radii (RS--in.) assumed for streamlines one through five at the 
stator exit plane, 

4, Same as (3,) but for the rotor (RR), 

5, Rotor loss coefficients, ZETAR, for incidence angle ratios of: 
0.0, 0.5, 1.0, -0.8, and -1.5, 

6. Rotor profile loss coefficients, ZETAPR, for incidence angle 
ratios of; 0,0, 0,5, 1,0, -0,8, and -1,5, 

7, Differences in stator gas outlet angles (DALF--radians) from 
the mean streamline angle, for streamlines one, three and five, 

8, Зате аз (7,) but for the rotor (DBET). 

9, Stator blade exit opening (AS--in,) for streamlines onc, three, 
and five, 

10, Same as (9,) but for the rotor (AR), 

ll, The number of sets of assumed conditions for which computations 
are to be made (NSETS). 

12, Assumed conditions for sets one through NSETS: stator mean as 
streamline exit Mach Number (AMS); inlet total pressure (PTO--psia); 
inlet total temperature (TTO--^R); rotor exit mean streamline Mach Number 
(AMR); wheel speed (RPM); and an indicator (NIND) which denotes whether 


assumed conditions except RPM are the same (NIND * O) or not (NIND = 0), 


176 


B, Outputs. 


Sample m is included at the end of the program, Since symbols 
used to label output ane Шалган” those used in the NS. printouts 
are self-explanatory for the most part. Printouts requiring clarification 
of meaning are given below in the order of the program output. 

l. І-( ) denotes the integrand of ( ). 

a XI = & 

3, PHI = È 

4, ARAT = a/an 

5, The SUMS 1-4, WSUM line gives the values of flow rate integrals 
between streamlines (in order) one and two, two and three, three and 
four, and four and five, followed by the sum of these values. 

6. STREAMLINE FLOW FRACTIONS, M = . Here, M denotes the stator 
(M = 1) or the rotor (M = 2). Numbers printed are the value of X = г/г„ 
followed by the flow fraction passing between the hub streamline and the 
streamline at X for streamlines two three and four in order. 

7. W-FRAC denotes flow fraction. 

8. The CONSTANT INTEGRAND 1-5 line of the ROTOR SOLUTION gives the 
values of the integrand for the angle part of Eq. (33) plus the entropy 
part 481 /4Х,, for streamlines one through five, 

9, INT2 is that part of the integrand of (33) for the remainder 
of the entropy tern. 

10. INT3 is that part of the integrand of (33) due to the equivalent 
enthalpy term. 

11. INT4 is that part of the integrand of (33) due to coriolis 
acceleration and rotor peripheral velocity. 

12. INT is the total integrand af (33)9 ores Ντ; re 


13, DSDXl-2 denotes dsj/dX,. 
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Units of printed values are: temperatures--°R; pressures--psia; 
angles--radians; velocities--ft/sec; radii--in.; and enthalpy terms-- 


BTU/1bm. 
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APPENDIX IV 
COMPUTER PROGRAM FOR EXHAUSTER ANALYSIS 
A, Inputs. 

Listed below are the inputs required for Program EJECTOR in the 
order of the READ statements. 

1. Тһе several values of polytropic efficiency assumed for the 
flow: through the turbine--ETAT; through the stator--ETAA; from the 
turbine hood to the ejector nozzle exit plane--ETAD; to the ejector 
nozzle throat--ETAN; from the ejector nozzle throat to the nozzle exit 
plane--ETAI; in the diverging portion (subsonic)--ETAC--and (supersonic) 
--ETAU; after a shock in the diverging portion--ETAL,  NSETS is the 
number of inlet conditions for which the analysis is to be done, 

2. Inlet conditions one through NSETS; total pressure (PTO--psia); 
and total temperature (TTO--°R). 

Β. Outputs. 

Sample output is included at the end of the program, Labels for 
the printouts are the symbols from the equations for the analysis. 

Units are: temperatures--^R; pressures--psia; velocities--ft/sec; and 


flow rates--lbm/sec. 
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APPENDIX V 
COMPUTER PROGRAMS PERTINENT TO FLOW RATE MEASUREMENTS 


A, Program ORIFICE. 

Inputs for this program are: the upstream pipe diameter (Dl--in.); 
the orifice diameter (D2--in.); the discharge coefficient, Ко, for the 
flange taps (AKF); the discharge coefficient, K*, for vena contracta 
taps (AKV); and the Reynolds Number for which the value of κ᾿ applies 
(REV). 

Outputs from the program are (see Table III, p. 102): the diameter 
ratio (BETA); the constant in the expression, X = C,w/z, (XC - C,); the 
constant in the expression «, for Y, (YC = Cy); the value of A for the 
flange taps (AF) and for the vena contracta taps (AVC); and the constant 
C for the flange taps (CFL) and for the vena contracta taps (CVC). 

B, Program FLOCAL. 

Inputs for this program, in the order of the READ statements are: 

l. The number of data points (N); atmospheric pressure (PA--in.Hg); 
reference manometer reading (TARE); and the base supply total pressure 
for the test (PUS--psia). 

2, For data points one through N: temperature (TN) and pressure 
(PN) upstream of the flow nozzle; differential pressure across the nozzle 
(HWN); temperature upstream of the orifice (TOR); pressures upstream of 
the orifice for flange taps (PFL) and for vena-contracta taps (PVC); and 
the differential pressures across the flange taps (HWFL) and the vena- 
contracta taps (HWVC). Temperatures are in °F, pressures are manometer 
readings (in.Hg) and differential pressures are in,H50. 

Sample output is included at the end of the program. Column head- 


ings are self-explanatory. 
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С. Program FLOCAL2. 


Inputs for this program in the order of the READ statements are: 

1. The total number of values of nozzle discharge coefficient to 
be read (N); the number of values of discharge coefficient to be read 
for test number one (Kl) and for test two (K2). 

2. Values of discharge coefficient one through N for flange taps 
(AK) and for vena-contracta taps (AK2). 

Sample output is included at the end of the program. The BEGIN and 
END values in the output specify which of the N values of discharge 
coefficient were used in the computations for the DATA GROUP. The 
column labeled POINT gives the number of the program's calculation 
subsequent to initial passage of the minimum error value of discharge 
coefficient, The desired value of discharge coefficient is obtained 
from the calculation which shows minimum error. 

D. Program LABRINT. 

No inputs are required for this program, Outputs are values of 
overall seal pressure ratio (BETA) and the corresponding values of the 
flow function (PHI = (0). 

E, Program LABLEAK, 

Inputs for this program are; 

1, The number of data points (N); ambient pressure (PA--in.Hg); 
and the number of throttlings of the labyrinth seal (Z). 

2. Test data for points one through N: reference manometer reading 
(TARE); temperature upstream of the square-edged orifice (T); pressures 
upstream of the orifice for flange taps (PFL) and the vena contracta 
taps (PVC); differential pressures across the orifice for flange taps 
(HWFL) and vena-contracta taps (HWVC)g turbine plenum temperature (TP); 


and turbine plenum total pressure (PPT). Temperatures are in μα 
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pressures are manometer readings (in.Hg) and differential pressures 
аге іп іп.Н20, 

Sample output is included at the end of the program, Column head- 
ings are self-explanatory. 
F. Program MINSQ2. 

Inputs for this program are; 

1, The number of sets of values of pressure and flow function to 
be read (N); and initial assumed values for the constants A and B, 

2. Values of pressure ratio (R) and flow function (W) for sets one 
through N, 

Sample output is included at the end of the program, Values of B 
and associated errors during final iterations are included in the output, 
Values of the flow function (м) as computed using the constant B, along 
with the experimental values, are printed for each test pressure ratio 
(R). In addition, values of the flow function computed from the para- 


bolic equation are given for pressure ratios from 0.15 to 0,95, 
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APPENDIX VI 
COMPUTER PROGRAM FOR TURBINE TEST DATA REDUCTION 


A. Inputs. 

Listed below are the inputs for Program TTTRPOS in the order of the 
READ statements, 

l. The number of data sets (NSETS). 

2. For each data set, the number of data points (NPTS), the run 
number (NRUN), the date of the run (NMOM) and (NDAY), ambient pressure 
(PA--in, Hg abs.), and ambient temperature (TA-- R). 

3. For each data point, one to NPTS, values of: the reference 
manometer reading (PTARE); the pressure upstream of the flow nozzle (PN); 
plenum pressure (PP); turbine inlet total pressure (PTO) and static 
pressure (PO); the stator exit plane pressures at the tip (PlT), the 
hub (PlH), and in the vicinity of the cover plate (PlI); and the turbine 
discharge pressure (P2). These pressures are manometer readings (in. Hg). 
If PlI has not been measured, the value of PlH is used for PII also. 

4, For each data point, one to NPTS, values of: flow nozzle 
differential pressure (HW--in.H50); the temperatures upstream of the 
flow nozzle (TN) and in the plenum (TP); the inlet total temperature (TTO); 
the axial tare force (AXTARE--1b); the axial force reading (AX) and the 
stator moment reading (AMOMS), both in kc.; rotor moment (AMOMR--in.lb.); 
and wheel speed (REV) in RPM, Temperatures are in °F, 

B, Outputs. 

Sample output from TTTRPGS is included at the end of the program. 
Column headings are the standard symbols for the turbine stations and 
are self explanatory. Temperatures are in °R, pressures are in psia.; 
angles are in degrees, and velocities are in ft/sec. The column labeled 


SPEED gives rotor mean streamline peripheral velocity, SE 
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